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Abstract 

We study a Lagrangian formalism that avoids double counting in effective field theories where 
distinct fields are used to describe different infrared momentum regions for the same particle. The 
formalism leads to extra subtractions in certain diagrams and to a new way of thinking about 
factorization of modes in quantum field theory. In non-relativistic field theories, the subtractions 
remove unphysical pinch singularities in box type diagrams, and give a derivation of the known 
pull-up mechanism between soft and ultrasoft fields which is required by the renormalization group 
evolution. In a field theory for energetic particles, the soft-collinear effective theory (SCET), the 
subtractions allow the theory to be defined with different infrared and ultraviolet regulators, remove 
double counting between soft, ultrasoft, and collinear modes, and give results which reproduce 
the infrared divergences of the full theory. Our analysis shows that convolution divergences in 
factorization formulas occur due to an overlap of momentum regions. We propose a method that 
avoids this double counting, which helps to resolve a long standing puzzle with singularities in 
collinear factorization in QCD. The analysis gives evidence for a factorization in rapidity space in 
exclusive decays. 
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I. INTRODUCTION 



Many problems of interest in quantum field theory have several momentum scales, and are 
efficiently treated using effective field theory (EFT) methods. One constructs a sequence 
of effective field theories which focus on one scale at a time. This greatly simplifies the 
calculations, partly because new symmetries emerge, and partly because Feynman graphs 
in each effective theory are much simpler to evaluate than the multiscale integrals of the full 
theory. More recently, theoretical methods have been developed which allow one to analyze 
field theories with several small momentum scales which are coupled by the dynamics. In 
these theories, it becomes necessary to treat the coupled momentum scales simultaneously 
within a single effective theory, rather than sequentially in a series of several different effective 
theories. Examples of theories with coupled scales are the soft-collinear effective theory 
(SCET) for energetic particles [fj], @, |], @], an d any non-relativistic theory where the kinetic 
energy is a relevant operator, examples being non-relativistic QED (NRQED) || and non- 



We wish to discuss an issue in the separation of infrared (IR) regions which appears at 
first to be a technical subtlety, but turns out to have important physical ramifications. It 
results in a tiling theorem for IR modes in quantum field theory. In the examples we discuss, 
it has to do with the proper treatment of soft modes in NRQCD/NRQED and of collinear 
and soft modes in SCET. In NRQED, the photon field of the fundamental QED theory is 
replaced by two fields, describing soft and ultrasoft (usoft) photons with energies of order 
mv and mv 2 respectively, where v <C 1 is the typical fermion velocity in the non-relativistic 
bound state. The soft and usoft NRQED gauge fields are A^(x) and A fJ, (x), where p is a label 
momentum of order mv, and k, the Fourier transform of x, is of order mv 2 H. The two fields 
describe photons with momenta p + k and k, respectively. In the special case that p = (the 
zero-bin), the soft photon becomes ultrasoft, and there is a double-counting of modes. To 
avoid double counting, the soft sector of the theory must have the additional constraint that 
p 7^ 0. This paper explores the consequences of implementing zero- bin constraints for soft 
modes in NRQED/NRQCD, and the analog for collinear modes in SCET. In loop graphs, for 
example, the sum over soft intermediate states should be ra f ner than the conventional 

J2 P - The difference, as we discuss in detail, is that conventional results have to be modified 
by zero-bin subtractions. 

The zero-bin subtraction solves a number of problems in NRQCD and SCET, and also 
resolves the long standing puzzle of divergent convolutions in QCD factorization formulas. 
We discuss several applications: 

1. Soft box graphs in NRQCD have unphysical pinch singularities in the energy integral, 



relativistic QCD (NRQCD) || % §, |, |I0|. 




(1) 
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which make them ill-defined, even in dimensional regularization. In previous computa- 
tions, it has been argued that these pinch singularities should be dropped in evaluating 
box graphs at any order in v |^, [11], |I2|| . Pinch singularities are also a problem for the 
method of regions A direct application of the method of regions for d A k leads to 



ill-defined integrals, so it was defined to apply to NRQCD only after first doing the 
energy integrals. The zero-bin subtraction modifies the soft box graphs so that pinch 
singularities are absent, and the graphs are well defined. 

2. The zero-bin subtraction automatically implements the previously studied pullup 
mechanism in NRQCD [[14], |15|1 , which was shown to be a necessary part of the defini- 
tion of this type of theory with multiple overlapping low energy modes. Through the 
pullup, infrared (IR) divergences in soft diagrams are converted to ultraviolet (UV) 
divergences and contribute to anomalous dimensions. 

3. There is a similar pullup mechanism at work in SCET for collinear diagrams. The 
anomalous dimensions of the SCET currents for endpoint B — > X s ^ and B — > X u lv 
were computed in Ref . |], [| from the 1/e and 1/e 2 terms. Some of these terms in the 
collinear graphs are actually infrared divergences. The zero-bin subtraction converts 
these infrared divergences to ultraviolet divergences so that IR-logs in QCD can be 
resummed as UV-logs in the effective theory. This formally justifies the results used 
for anomalous dimensions in these computations, and in subsequent work for other 
processes with similar anomalous dimensions eg. J[£| [I7[ [TJ, [L9|, |2y, |2Ifl . 



4. In high energy inclusive production such as 7* — > qqg, there is a potential double 
counting at the edges of the Dalitz plot in SCET, which is resolved by properly tak- 
ing into account the zero-bin in both fully differential and partially integrated cross 
sections. 

5. As a by-product of our analysis we give definitions for NRQCD and SCET that are 
independent of the UV and IR regulators. We also demonstrate a link between power 
counting and reproducing infrared divergences in the EFT. Exploiting this link, we 
demonstrate that a choice exists for the degrees of freedom in NRQCD and SCET 
which is complete, covering all infrared regions for a broad class of physical situations. 
For these cases no new modes are required at any order in a s or in the power expansion. 

6. In high energy exclusive production, such as 7* — > up or 7* — > hit, there are un- 
physical singularities in convolution integrals of some hard kernels with the light-cone 
wavefunctions (j) n (x). For example 

"d.tM (2) 
x 
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which is divergent at x — > if ^(x) vanishes linearly as x — > 0. The same is true 
for exclusive light meson form factors at large Q 2 , as well as processes like B — > 7r£z/ 
and i? — > 7T7T for ^ ^> Aqcd- The zero-bin subtraction removes the singularity, 
and induces a corresponding UV divergence. After renormalization we have a finite 
convolution with (j) n (x), and the kernel behaves as a distribution we call 0: 

\ x *M^[\ x p^ = [\ x mzM^m +k<ao . (3) 

x Jo \ x )0 JO x 

The (^Tr-type term involves a In (£!„•) and is induced by UV renormalization in rapidity 
space. These terms are discussed in the body of the paper. Thus using SCET finite 
amplitudes are obtained for apparently singular hard scattering kernels. 

7. The zero-bin procedure gives insight into factorization formulas which separate modes 
in rapidity space rather than by scale separation in their invariant mass. For example, 
two hadrons both built of non-perturbative modes with p 2 ~ Aq CD can have their 
modes factorize by being in different corners of phase space or rapidity. We discuss 
this by applying the zero-bin technique to the formulation of degrees of freedom in 
SCET n with Wilson lines on the light-cone. The separation between soft and collinear 
regions is controlled by perturbation theory with dependence on a rapidity parameter. 

Examples of processes where endpoint singularities in convolution integrals have been en- 
countered include the pion form factor at large Q 2 and subleading twist | 22|| , the p — it 



form factor [23], the B — > ittv form factor form factor terms in B — > txtx [25], the 



Pauli nucleon form factor F2 [26], color-suppressed 5-decays involving light isodoublets, 
such as -B — > D°K° and B° — > D S K~ ||27||, and annihilation contributions in two-body 



£>-decays [^8], ^9], |30|. Endpoint singularities also appear in non-exclusive processes such 
as semi- inclusive deep inelastic scattering at low transverse momentum [[H], [32]. When the 
zero-bin procedure is applied to these cases, individual pieces of the amplitudes exhibit 
dependence on a rapidity parameter. 

In the work of Collins and Soper |3B], factorization formulas involving a rapidity 



parameter were derived for fragmentation in e + e~ — > A + B + X where A and B are 
hadrons. The separation of degrees of freedom in SCETn gives finite amplitudes that appear 
to indicate that rapidity dependent factorization also occurs in other two-hadron processes, 
including purely exclusive ones. Our finite amplitudes are shown to be a direct consequence 
of defining the degrees of freedom in SCETn carefully. We will not give a complete derivation 
of an exclusive rapidity factorization formula here, because in our analysis we will make the 
simplifying assumption that the renormalization of rapidity space effects and invariant mass 
effects can be carried out independently. 

It is useful to have a physical understanding of why resolving the double counting issue 
also resolves the singularity problems. A hint comes from the fact that neither the pinch 
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singularities or endpoint singularities are present in the full QCD computations, as pointed 
out for SCET in |35| in the context of B — > iriu. For case 1) above, the soft pinch singularities 
are removed by the kinetic energy of the quarks, l/[k° — k 2 /(2m)]. In case 7) above, the 
endpoint singularities are soft limits of the full theory diagrams and are removed by Aqcd; 
or in non-relativistic systems by the binding energy. An improper interpretation of the 
singularities can occur by not being careful about taking double limits. In 1) the singularity 
arises from first taking k° ^> k 2 /(2m) and then k° — > 0. The k° k 2 /(2m) limit gives 
soft quarks which can not properly describe the potential-quark region where k° ~ k 2 / (2m). 
Likewise, in 7) the endpoint singularity comes from first taking k~ 3> k±, k + and then taking 
k~ —>■ 0. The collinear particles obtained from k~ 3> k±, k + do not properly describe the 
soft particle region where k~ ~ k±,k + . To avoid double counting we must ensure that 
the soft quarks do not double count the potential region in case 1), and that the collinear 
quarks do not double count the soft region in 7). In the effective theory implementation 
of the expansion of QCD these singular limits are properly described by other degrees of 
freedom. Once we avoid the double counting, the unphysical singularities never appear 
because the potential limit of the soft quarks and the soft limit of the collinear quarks are 
rendered harmless. It should be emphasized that the pinch and endpoint singularities we are 
discussing are unphysical artifacts of certain approximations, and are reflected by double 
counting problem that must be fixed in the Effective Field Theory (EFT). They are not 
the same as the classification of physical IR divergences from the Landau Equations (see 
eg. |J7], [SB]) which go by similar names. The true IR structure of the full theory is properly 
reproduced by contributions from the full set of EFT degrees of freedom. 

It is important to note that the zero-bin subtractions avoids double counting independent 
of the choice of UV and IR regulators in effective theory computations. Our implementation 
of zero-bin subtractions is unique up to possible finite scheme dependent contributions, and 
provides an explicit connection to methods which introduce hard factorization cutoffs. It 
also provides a definition of the modes in these effective theories independent of dimensional 
regularization (and with some work could be used for example to take the cutoff formulation 
of SCET described below and implement it on the lattice). Since the proper formulation of 
an EFT should not depend on the choice of IR regulator used in perturbative computations, 
this is not surprising. 1 Physical results in QCD are IR finite with divergences removed by 
quantities like Aqcd ; binding energies, or cancellations between real and virtual diagrams, 
and the same is true in the EFT. 

The outline of the paper is as follows. In section |I| we give a brief introduction to 
NRQCD and SCET which serve as our main examples. The tiling of IR regions with modes 



The choice of IR or UV regulators can make it more difficult to perform the power counting, for example 
by leading to integrals that are not homogeneous and require power counting violating counterterms to 
give back a power counting for the renormalized graphs. 
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is discussed in section pTTj , and the zero-bin subtraction is formulated in section [IV]. Our 
discussion is in the context of NRQCD and SCET, but is general enough to be readily 
adapted to other physical situations. In section |V| we give examples in NRQCD to show 
that the zero-bin removes pinch singularities, implies the pullup mechanism, and avoids 
double counting problems. In section |VI| we give examples in SCETi which demonstrate the 
regulator independence of the zero-bin method, and the removal of double counting in loop 
integrals and in inclusive phase space computations. Finally, in section |V1I| we give examples 
in SCETn which is formulated with zero-bin subtractions and only soft and collinear modes. 
The zero-bin subtraction resolves the endpoint singularity issue in exclusive processes to 
leave finite amplitudes, and require the introduction of a rapidity parameter. Conclusions 
are given in section |VIII| . 



II. NRQCD AND SCET 

NRQCD is an effective theory for non-relativistic quark-antiquark (QQ) bound states, 
where the typical relative fermion velocity, v, is small, ?; < 1. The relevant scales in NRQCD 
are the quark mass, m, momentum p ~ mv , and energy E ~ mv 2 , with E <C p <C m. The 
energy and momentum are not independent; they are coupled via the quark equation of 
motion, 2mE = p 2 . 

SCET describes the interaction of energetic particles; examples include the inclusive decay 
B — > X s 7 at large E 1 via the partonic decay b —>■ s'j, inclusive jet production, or exclusive 
semileptonic decays such as B — > nlv. It is convenient to orient the coordinate system in the 
direction of the energetic jet or hadron, for example by introducing null vectors n = (1,0,0, 1) 
and n = (1, 0, 0, —1), and use light-cone coordinates with p + = n ■ p, p~ = n- p for any four- 
vector p. Energetic particles moving near the n direction have momenta p~ ~ Q, p + ~ QX 2 
and p± ~ QX, where A 1 and Q S> Aq CD . Q is the large energy scale, and is of order m h 
for £?-decays. Often the choice A ~ a/Aq CD /Q is made for inclusive processes in SCET I; but 
parametrically larger choices for this small parameter are allowed. For exclusive processes 
in SCETn we use r] rather than A for the expansion parameter to avoid confusion, since here 
r] ~ Aqcd/Q- The p + , p~ and p_i_ scales are coupled by the on-shell condition p + p~ = p\, 
and despite the hierarchy p~ ^> p±_ ^> p + , the effective theory must simultaneously deal 
with the low energy scales associated with p± and p + . 



A. Comparison with the NRQCD Method of Regions 



One can evaluate Feynman integrals in NRQCD using the method of regions [|13|] (also 
called the threshold expansion), which divides up an integral into hard, soft, potential and 
usoft contributions based on scaling for the loop momentum. The idea is that in dimensional 
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regularization the sum of these contributions exactly reproduces the full theory diagram, so 

A iu \ Pi ) = n j d % ^ = e n \ d % F(e) ^ k i) ■ ( 4 ) 

j regions £ j 

Dimensional regularization is required here because there are cancellations between UV 
and IR divergences from different regions, which occur if em = euv = e- Eq. (f|) does 
not define an effective field theory, but it is sometimes taken as a way of defining EFT 
contributions to amplitudes, by demanding that in dimensional regularization each mode in 
the effective Lagrangian should reproduce a term from a region on the RHS. Although Eq. 
is quite powerful, a few points must be treated carefully: i) The division of regions is gauge 
dependent. 2 ii) The requirement that scaleless integrals be set to zero, 1/euv — 1/eiR = 0, 
does not allow all UV divergences to be treated by counterterms, nor a verification that 
every IR divergence has a correspondence with the full theory, iii) In Eq. (^) one must sum 
over all possible momentum routings in loops to determine the relevant regions (or consider 
the scaling of all combinations of loop momenta and external momenta). This is because it 
is individual propagators in the EFT that belong to a region rather than the loop momenta. 

In the remainder of this section we will explore how the terms (1/euv ~~ l/ e ni) allow a 
residual freedom in associating amplitudes with degrees of freedom beyond that in Eq. (f|). 
The treatment of these terms effects the correspondence of the EFT modes with physical 
regions. Lets consider a one-loop graph in NRQCD with contributions from different regions. 
The soft contribution depends on external soft scales such as the momentum transfer r, and 
has the (schematic) form 

A B 

/soft = + — + /(r,//). (5) 

As the momentum in the soft graphs vanishes, the graph matches on to an usoft diagram, 
with the (schematic) structure 

B C 

^usoft = 1 \-g(E,/j), (6) 

euv eiR 

where the coefficient B is the same as in Eq. @ . The usoft graph depends on external usoft 
scales such as the energy E. The IR divergences in the usoft sector are true IR divergences. 
They arise if one is computing an IR divergent quantity such as an on-shell Green's function, 
but cancel in measurable quantities such as physical scattering cross-sections and bound state 
energies. 

The IR divergence in the soft graphs and the ultraviolet divergence in the usoft graph 
are at the intermediate scale mv and cancel each other; they are not true divergences of the 



2 A well known example in NRQCD is applying the method of regions to potential and soft contributions. 
Another example is the division between (E, p) ~ (m, to) and a new region (E, p) ~ (m, mv) that shows 
up at fourth order in the v expansion in Coulomb gauge. 
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FIG. 1: Comparison of two setups for the soft and usoft contributions to an NRQCD Feynman 
graph. In a) the 1/e divergences at the intermediate scale mv cancel between the soft and usoft 
contributions. In b) there are no IR divergences at this intermediate scale. 



theory. Since in the method of regions a rule is applied that scaleless integrals are set to 
zero, one is free to consider the B/e terms canceling in the sum of Eqs. (flU) to give 

A C 

/soft + /usoft = 1 H /(r, fi) + g{E, fj). (7) 

This interpretation leads to the picture shown in Fig. |l]a, where the IR soft effects and 
UV ultrasoft effects meet at mv. Thus this particular mapping of the method of regions 
with effective theory amplitudes does not lead to simultaneously having degrees of freedom 
for soft and usoft gluons at a scale /i, but instead in the renormalization group the soft 
contribution covers the momentum region between m and mv and the usoft contribution 
covers the region between mv and mv 2 . 

In an effective field theory, the low energy effective Lagrangian is usually required to 
reproduce all the IR effects of the original theory from the start. In NRQCD, this plus 
a manifest power counting in v requires that the effective theory include soft and usoft 
degrees of freedom at the same time to reproduce the dependence of the full theory on the 
momentum transfer r ~ mv and the energy E ~ mv 2 , which are considered IR scales at 
H — m where the effective Lagrangian is constructed, as found in Ref. [[|. This happens 
because the power counting links the scales mv and mv 2 . A similar result holds for SCET — 
the SCET Lagrangian at the hard scale Q must simultaneously include collinear and usoft 
degrees of freedom to correctly reproduce the IR behavior of QCD with a power counting 
in A [§], ID- In the effective theory, counterterms must be added for all 1/euv divergences, 



including those from scaleless integrals as is familiar from the study of HQET j3£], [40], |41 
The counterterm structure can not depend on the choice of IR regulator. 

The above line of reasoning gives the picture Fig. [Tft> , rather than Fig. |l]a. In NRQCD 
this picture has been implemented in the past by a procedure referred to as a pullup mech- 
]l4l , [l5| . One modifies Fig. |l]a by pulling up the usoft modes to the scale m, and then 



amsm 



subtracting the usoft contribution from the soft contribution in the region between m and 
mv. The picture now looks like Fig. |l]b. This modification leaves the usoft contribution 
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Eq. (|6|) with the same form, but it must now be included in the theory between the scales 
m and mv 2 . The usoft integral in the region between m and mv has the form 

B B 

(8) 



so it corresponds to a (1/euv — V e nt) that we are free to move between amplitudes from 
the point of view of Eq. (§]). When Eq. (|]) is subtracted from Eq. (H) it gives the modified 
soft contribution 

/soft = — + — + /(r,/x). (9) 

The IR divergent parts of the soft graphs have been removed. There is no need for any 
cancellation of divergences between the soft and usoft degrees of freedom, and they can be 
renormalized separately. This is important, because it is known that the soft and usoft 
degrees of freedom have independent coupling constants, and the cancellation of divergences 
between the two can be problematic in renormalization group improved perturbation theory. 
One of the consequences of the pullup mechanism is that the anomalous dimensions in the 
soft sector are given by Eq. (|9]), and are proportional to A + B. As a short cut one can 
compute the anomalous dimension from the original form in Eqs. by treating the IR 

divergence as though it were a UV divergence in Eq. (|5[), and taking the UV divergence in 
Eq. ([]) to be at the hard scale. This is the procedure that has been followed in previous 
NRQCD and NRQED computations g 0, [Lj, |TJ, g^, (£| g§ g5|, and is known to be nec- 



essary to correctly reproduce the high-order logarithmic terms in Lamb shifts and hyperfine 
splittings which are determined from independent fixed order QED computations. We will 
see that the zero-bin subtraction automatically gives the final result implied by the pullup 
mechanism, so that the soft integrals with the zero-bin subtraction have the form in Eq. @, 
rather than that in Eq. (|5|). Thus the zero-bin formulation no longer requires implementing 
a pullup by hand. 

An alternative setup has been explored in Ref. [46 , 47 , 48|, 49| , which keeps Fig. p]a 



but postpones obtaining power counting in v. In this setup one starts by matching onto 
a purely soft theory with a power counting in 1/m and constrains states to E = so the 
IR divergences are properly reproduced. The power counting in v is not yet manifest in 
this theory. The soft theory is then matched at a soft cutoff scale Q s onto an usoft theory, 
pNRQCD, where E ^ and the velocity power counting is restored. The key is to maintain 
VL S as a variable in the usoft theory so that when required it can be run down from the scale 
m including the required correlation with the usoft scale. While this setup seems to differ 
from the pullup, in several cases it has been shown that both methods give equivalent final 
answers for observables with log summation. 

Note that so far an analogous setup that would avoid simultaneously introducing usoft 
and collinear modes does not exist for SCETj. It is known that one can avoid simultaneously 
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introducing hard-collinear p 2 ~ QAqcd and collinear p 2 ~ Aq CD modes, by matching SCETi 
onto SCET n as discussed in Ref. p5f . This postpones obtaining the final power counting 
until one matches onto SCETn, and the (hard-collinear)-(collinear) setup is similar to the 
soft-usoft setup in pNRQCD. Just as in pNRQCD one must in general maintain the matching 
scale as a free parameter with this method of matching SCETi onto SCETn- Alternatively, 
one can match QCD directly onto SCETn @, |50 . 



B. Summing Logarithms 

The dynamical relations p 2 = 2mE and p + p~ = p 2 ± have implications for the summation 
of logarithms using renormalization group evolution. In NRQCD and NRQED, one must 
simultaneously run from m — > p and m — > E in the soft and usoft sectors of the theory ||. 
This is implemented in practice by using the velocity renormalization group ||. Graphs 
in the theory are evaluated using two different fi parameters, fis for soft and potential 
loops, and nu for usoft loops. One then sets [is = Tnu, jiu = mu 2 , and runs from v = 1 



to v = v. This procedure is also referred to as one-stage or correlated running [51], and 
corresponds to Fig. |l]b. In NRQED, this correlated running is required to correctly compute 
the a 8 In 3 a Lamb shift and a 7 In 2 a hyperfine splittings for Hydrogen and positronium, as 
well as the a 3 In 2 a positronium widths f43| , [5lH . It has also been shown to be necessary to 
properly implement counterterms in sub divergences at three-loops fP], |52fl . In some cases, 
correlated running is not essential, and one can follow an alternative procedure called two- 
stage or uncorrelated running, in which fi$ is scaled from m to mv, and flu from mv to mv 2 , 
corresponding to Fig. |l]a. 3 The summation of logarithms for the 1/m 2 QCD potentials can 



be done with or without the correlated running j|2], [47], both methods give the same 
result. 

So far no examples where correlated running is essential have been encountered in SCET. 4 
This does not mean that correlated running is not necessary in SCET. In NRQED, correlated 
running is first required for computations of recoil corrections at order v 3 (m e /m p terms in 
the Lamb shift), because it is at this order that the potential and usoft divergences are tied 
together. Until this order, both correlated and uncorrelated running give the same result, 
and it is the Lamb shift computation which shows that, in general, one should use correlated 
running. Thus, it is very likely that only the running of subleading factorization formulas 
in SCET will demonstrate the manner in which correlated running is required. 



3 See Ref. |5lJ for the precise relation between these two methods. The anomalous dimensions have different 
definitions in the two approaches, so the single In terms agree. There is a difference only for In 2 and higher 
terms if the anomalous dimension does not factor into separate soft/potential and usoft pieces. 

4 See the discussion in Ref. Q on the equivalence of correlated and uncorrelated running at leading order 
for B — > X s -f and e + e~ — > J/tpX. 
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III. INFRARED MODES IN NRQCD AND SCET 



In its region of validity, an effective field theory needs to systematically reproduce the 
IR structure of the full theory order by order in its power expansion. In both NRQCD and 
SCET, a strict interpretation of this requirement makes it necessary to include distinct fields 
for different moment regions of the same physical particle. These distinct fields have different 
power counting. Multiple gluon fields were first introduced in Ref. |54| for potential and 
usoft gluons in NRQCD, with energy and momentum scaling of order (E ~ mv 2 ,p ~ mv) 
and (E ~ mv 2 ,p ~ mv 2 ), respectively. We will not introduce fields for potential gluons since 
they are not propagating degrees of freedom. In NRQCD it is also necessary to introduce 
soft gluons with momentum scaling (E ~ mv,p ~ mv) JIB] , |55fl . We use the NRQCD 
Lagrangians defined as in |5|, |M| with potential quarks and usoft gluons/quarks, and also 



soft quarks. We treat soft gluon vertices with a soft-HQET effective Lagrangian, rather than 
integrating out the soft quarks as in Ref. |j4]. At two-loops and beyond it is important to 



keep track of the i0 + in the soft quark propagators l/(v ■ k + i0 + ) and leaving them in an 
action facilitates this. In SCET one requires both collinear and usoft gluons in a theory often 
called SCETi or collinear and soft gluons in a theory called SCETn • We use the Lagrangians 
from Ref. |4j] for these theories. 

One might expect that introducing multiple fields for the same particle would lead to 
double-counting problems. In constructing effective theories, one needs to know not only the 
power counting for the degrees of freedom, but also understand the range of scales for which 
these modes are included in the effective Lagrangian. If degrees of freedom overlap in some 
region of momentum space, an understanding of how their definitions avoid double counting 
is necessary. In the effective theories we study, loop integrals are dominated by external 
momenta by construction, so the power counting guarantees that fields give contributions 
that can overlap only in UV or IR limits. 5 As a simple toy example consider Fig. |^. We 
imagine that there are two relevant momentum coordinates pi and P2, and that physically 
there are four interesting sets of momenta labeled q a , qb, q c , and qd which could be set 
kinematically or by bound state dynamics. The hard cutoffs Ai and A 2 distinguish the 
momentum regions dominated by these g's. 

Consider first the simplified case where we ignore the p<i axis, and only have q a and 
This situation applies to many physical problems, including that of integrating out massive 
particles like the W^-boson or fe-quark. Here denotes hard fluctuations that are integrated 
out into Wilson coefficients C, while q a denotes low energy IR modes. To simplify the 
renormalization group evolution and leave symmetries unbroken, it is convenient to trade 
Ax for a scaleless regulator such as dimensional regularization. The low energy theory has 



5 The statement actually holds for any renormalization procedure that respects the power counting, or at 
worst requires power counting violating counterterms. 
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FIG. 2: Toy model to illustrate the scales captured by an effective theory with multiple low energy 
modes (here q a , qb, and q c ). 



C = C(qb,fi) where fi is the dimensional regularization parameter. The effective theory 
for q a with a scaleless regulator takes A x — ► oo, and thus overlaps with the hard region of 
momentum space, but only through ultraviolet effects. In the effective theory this double 
counting is removed by UV counterterms oc 1/e, as well as through finite terms in the Wilson 
coefficients. 

Now consider both the p\ and P2 axis. Here qd denotes hard fluctuations, and q a , <#,, and q c 
are all low energy modes because they border a region where one or both of px,2 can be zero. 
Using dimensional regularization for each of these modes effectively takes A 12 — > oo for q a ; 
Ai — > 0, A 2 — > oo for 5b ; and A 2 — > 0, Ai — > oo for g c . Again there is double counting in the 
UV which is taken care of by counterterms. However, there is now also a double counting 
in the IR. For example, taking Ai — > for qb runs into the region for q a . The Ai — > limit 
is necessary for any scaleless regulator. In fact this limit is actually important physically, 
because we would like to define qb without reference to q a in order for it to be possible that 
contributions from q a and qb can factorize into independent well defined objects. This double 
counting in the IR is removed by the pullup mechanism, which as we will demonstrate, is 
a consequence of a proper treatment of the zero-momentum bin, namely p 2 = for q c and 
Pi — for qb. In pure dimensional regularization these zero-bin contributions amount to a 
correct interpretation of terms cx 1/euv — V e rR- F° r example, 1/euv terms must be canceled 
by counterterms while 1/erR terms match up with IR divergences from the full theory. 

We now turn to realistic effective theories and their zero-bin's. For NRQCD the modes 
are shown in Fig. |3|a. The hard scale is E ~ m or p ~ m. We have propagating soft and usoft 
gluons/light quarks with power counting E ~ p ~ mv and E ~ p ~ mv 2 respectively, and 
potential heavy quarks with E ~ mv 2 and p ~ mv, and soft heavy quarks with E ~ p ~ mv. 
Hard cutoffs A, A 1; and A 2 have been introduced to facilitate the discussion. Here A ~ m 
is an ultraviolet scale below which one uses the effective theory, and A i)2 divide up the low 
energy modes. These cutoffs will be removed exactly as in the toy example above. First 
consider the soft and usoft gluons. In a theory with both present, the energy /momentum 
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FIG. 3: a) Scales and momentum modes for nonrelativistic field theories like NRQCD. Here s, p, 
and u denote soft, potential, and usoft respectively. 



regions we want them to cover are 



soft gluons: 
usoft gluons: 



A > mv > tt s 
A»mt) 2 >n M , 



(10) 



where Q s and Q u are soft and usoft scales denoted by the box and star in Fig. || These 
scales are usually set by external variables such as the momentum transfer or energy of the 
quarks, or by the nonperturbative scale Aq C d- In order to reproduce all possible IR effects 
associated with E ~ mv 2 , the EFT necessarily must have usoft modes just below the UV 
scale A. The division in Eq. (|T0|) also implies that the UV divergences associated with soft 
and usoft modes contribute to the same anomalous dimension. This setup corresponds to 
the result obtained from the same limits for the cutoffs A 12 as discussed in the toy model 
above. Now double counting must be avoided in the region between A and Q s , between the 
IR of the soft gluons and the ultraviolet of the usoft gluons (i.e. when virtual momenta for 
the soft gluons becomes comparable to virtual momenta for the usoft gluons in loops). 

A line of reasoning that gets close to seeing how this is achieved is to start with soft 
gluons over the interval A ^> mv ~ Ai and usoft gluons over Ai ^> mv 2 ~ Q u , where we 
set A 2 = Ai- Here Ai is considered as an intermediate factorization scale ^7], fjg| . Next 
one adds a pullup contribution to the usoft gluons and simultaneously subtracts the same 
contribution from the soft gluons as described in Refs. [|14], [15| . This pulls the upper limit 
Ai for the usoft gluons all the way up to the scale A, while at the same time avoiding double 
counting between m and mv. The pre-pullup setup is shown in Fig. |l]a and post-pullup in 
Fig. |l]b. Taking into account the zero-bin gives the post-pullup setup. With the pullup, 
the contributions from soft and usoft gluons are now as in Eq. fllPf) . If the cutoffs had been 
swapped for a scaleless regulator, this would amount to a proper interpretation of 1/e poles 



as discussed in Ref. [15 
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In this paper we show that Eq. (|10D is obtained automatically by carefully considering the 
zero-bin in the NRQCD effective Lagrangian. The usoft gluons are defined all the way up 
to A, and the subtraction for the soft gluons is associated with properly removing their zero 
momentum bin. Thus the effective Lagrangian gives the full procedure for the evaluation 
of Feynman graphs and is decoupled from the additional choice of which regulators to use. 
With the zero-bin taken into account there is no longer a need to implement a separate 
pullup. 

A second type of division between potential and soft IR modes in NRQCD is also shown 
in Fig. |3|a. In this case the distinction is solely in the energy variable p°. Potential gluons 
are not propagating degrees of freedom since they have p° ^ p ; and so potential gluon fields 
should not be introduced (they would have problems with gauge invariance for example). 
The matching of soft gluons onto four-quark operators with potential coefficients can be 
thought of in a similar manner to integrating out a massive particle |56| . On the other 
hand a zero-bin subtraction is necessary to distinguish soft and potential quarks. Their 
momentum regions and propagators are 

soft quarks: (p° : A > p°~?m> > Q s ) ; propagator: — , (11) 

[p u + i0 + \ 

potential quarks: (p° : A > p° ~mi? 2 > Q u ) ; propagator: ^ • 

Without the zero-bin, double counting occurs when p° — > in the soft propagator, and this 
reveals itself through the presence of pinch singularities in soft loop diagrams with quarks 
and antiquarks, which have the form 

do _ (12) 

P b° + iO+][-p° + iO+]...' 1 ; 

Here we will show that the zero-bin subtraction removes all pinch singularities in the soft 
regime. At the same time it avoids double counting of soft and potential contributions. It 
is useful to recall that the same is not true in the method of regions ||13|| , where the method 
used to avoid the pinch singularities is to consider the expansion only after doing the p° 
integral by contours. 

Note that in NRQCD there is no issue of a possible potential-usoft overlap for quarks 
since propagating usoft quarks are light quarks and therefore of a different flavor from the 
potential quarks. 

So far the discussion was for NRQCD, but similar logic holds for SCETj. In this case 
the hard scale is set by Q and the expansion parameter is A. Below the scale Q we have 
usoft and collinear gluons with power counting (p + ,p~ ,p x ) ~ Q(X 2 , A 2 , A 2 ) and Q(X 2 , 1, A) 
respectively. These gluons cover the regions 

usoft gluons: (p" : A > QX 2 > Q u ) ; (p ± : A > QX 2 > Q u ) , (13) 
collinear gluons: (p~ : A > QX° > Ql) ; (p L : A >QA > Q° ± ) , 
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Qrf Ai Gri A 2 er f A P+ 

FIG. 4: a) Scales and momentum modes for SCETj. Here cn, cn, and u denote collinear-n, 
collinear-n, and usoft modes respectively and A ~ \/Aqcd/Q- b) Scales and momentum modes 
for SCETn- Here cn, cn, and s denote collinear-n, collinear-n, and soft modes respectively and 
r} ~ Aqcd/Q- 



and a common region for p + \ A ^> QA 2 ^ Q u . These regions are denoted by a circle and star 
in Fig. |]a, where we have also included a second collinear region cn for later convenience. 



In Eq. (|13j) we have used a common UV scale A and a common usoft IR scale f2" in the — 
and _L components which have power counting ~ A 2 . In SCETi we see that we must avoid 
double counting in the region A to Q c _ for the p~ momenta, and the region A to Q c ± for the 
p 1 - momenta. 6 As in NRQCD, this is achieved by a proper treatment of the zero momentum 
bin for collinear fields (which implements a pullup in SCET). 

Finally, we can consider the theory SCETn which has soft and collinear gluons with power 
counting p~ ,p L ) ~ QiViViV) an d Q(v 2 AiV) respectively, for r\ ~ Aqcd/Q- Here the 
regions are 



soft gluons: 
collinear gluons: 



(p- : A > Q V > VI s ) ; (p A 
(p- : A > Q V ° > ; (p H 



A > Q V > VI s ) , 
A > Qrj 2 > n%) , 



(14) 



and for convenience a common region for p- 1 : A ^> Qr] > Q s . Both soft and collinear 
modes describe non-perturbative fluctuations close to the mass shell, p 2 ~ Q 2 1] 2 ~ Aq CD . 
Interactions between these modes are offshell by an amount, p\ c ~ Q 2 r], and dependence 
on this momentum is integrated out, appearing in the coefficient functions for mixed soft- 
collinear operators. Here double counting in SCETn occurs when a collinear momentum 
overlaps the soft region, and when a soft momentum overlaps the collinear region. This case 
differs from our discussion of NRQCD and SCETj because here the overlapping modes have 



For typographical convenience we use f2_ although a superscript, , would be more appropriate. 
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the same p 2 , but differ in their rapidity y, or more conveniently their value of 



C = ^ = P ¥ • (15) 

The n-collinear modes have ( p 3> 1, the soft modes have ( p ~ 1, and the n-collinear modes 
have ( p <Cl. Consider a process for which n-collinear modes are irrelevant. Double counting 
is avoided by a proper treatment of the zero-bins: the = bin" for n-collinear modes 
and the u p + = bin" for soft modes. However, here double counting of a physical IR region 
in QCD requires a correlated change in the + and — momenta: p~ gets small while p + 
gets big for collinear, and p~ gets big while p + gets small for soft. The implementation of 
zero- bin's in SCETn is discussed further in section IV11. 



The treatment of the zero-bin ensure that the double counting is removed in the infrared, 
and that the overlap in the ultraviolet is properly handled by renormalization in the effective 
theory irrespective of the choice of regulator. If dimensional regularization is used to regulate 
both the IR and UV, then the added contributions are scaleless loop integrals that appear 
to be zero. In logarithmically divergent integrals, this occurs because the integral is the 
difference of UV and IR divergences, = 1/euv — l/ e iR, since there is only one e. If 
these added contributions are ignored, then one must be careful to properly interpret the 
divergences as UV or IR. This conversion of IR to UV divergences has been used implicitly 
in much of the NRQCD and SCET literature. However, if one wants to fully understand the 
physical significance of certain divergences or use another regulator in the UV or IR then 
explicitly including the subtractions discussed here is necessary. 



IV. ZERO-BIN SUBTRACTIONS 



We start by reviewing how the relevant momentum scales are separated in the effective 
theory using labeled fields. In NRQCD, one first removes the large mass m of the quark (and 
antiquark) from the problem just as in heavy quark effective theory (HQET) p7| . The total 
momentum of the quark, P M , is written as the sum P M = mv M + g M , where v M = (1,0,0,0). 
This subtracts m from all the energies. The residual momentum g M is much smaller than m, 
and contains the non-relativistic energy E ~ mv 2 and momentum p ~ mv of the particle, 
where v is a scaling parameter of order the typical relative velocity between the heavy quarks. 
This mixes different powers of v. As shown in Ref. |§, it is useful to make a further division 
of g M , g M = p^ + k^ 1 , where p M is of order mv, and k^ is of order mv 2 . This second separation 
allows the power counting in v to be manifest in the effective theory. The break-up of g M is 
shown schematically in Fig. ^a. One breaks momentum space into a discrete variable p of 
order mv, and a continuous variable k^ of order mv 2 . Often k is referred to as the residual 
momentum. The discrete label p does not have a time-component for quarks, since the 
energy is of order mv 2 . The entire g M momentum space is covered by integrating over k^ 
and summing over the labels p. Quarks are described by fields ip p (x), with an explicit label 
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FIG. 5: Label and Residual momenta for a) NRQCD quarks and b) SCETj. In both cases p 
denotes a large momentum, and labels a particular box, whereas k is a small momentum, and gives 
the final momentum location relative to the reference momentum point in the box labeled by p. 



p, and the momentum is the Fourier transform of x. Similarly, soft gluons with energy 
and momentum of order mv are described by gauge fields A£(x) with a four- vector label 
p^. Ultrasoft gluons with energy and momentum of order mv 2 are described by gauge fields 
A^ix). Those unfamiliar with how the field theory with label and residual momenta works 
are referred to Ref. pi] or the example in the next section. 

An analogous procedure was applied to SCET in Ref. @, ||]. In SCET, one breaks up 
the collinear momentum into a label and residual momentum, P M = p^ + k^. The label p 
contains the Q and QX pieces of the momentum, and the Fourier transform of the coordinate 
x is the QX 2 part of the momentum k. Unlike NRQCD, the collinear interactions can still 
change the large label momentum ~ Q. In a theory referred to as SCETi, one has collinear 
quark and gluon fields £ n ,p(x) an d A% (x), where the label p has minus and _L components, 
and usoft quark fields A^(x) which have energy and momentum of order QX 2 . The SCETj 
decomposition is shown schematically in Fig. [|b. In a theory called SCETn one has collinear 
quarks and gluons, £„ iP , A£ p with p~ label momenta, and soft quarks and gluons, q s ^, A^ k 
with k + label momenta. Here soft fields have plus-momenta much bigger than their collinear 
counterparts, and collinear fields have minus-momenta larger than their soft counterparts. 

The quark and soft gluon fields in NRQCD and the collinear quark and gluon fields in 
SCET will be referred to as labeled fields. As is clear from Fig. [|, labeled fields must have 
a non-zero value for their label, i.e. they must be outside the zero-bin. Otherwise, they 
cover the same momentum region as the usoft fields. This is implemented at the level of the 
effective theory Lagrangian by requiring that all labeled fields in the Lagrangian have a non- 
zero value for the label. Terms in the effective theory Lagrangian have sums over the field 
labels, and all such sums are over non-zero values of the labels. One then has to carefully 
derive the rules for effective theory graphs including this constraint on the labels. 

In loop integrals, one finds expressions which involve a sum over labels and an integral over 
the residual momentum, which can be converted to an integral over the entire momentum 
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space [§], 

[dk—* [dp, (16) 



p 



as is clear from Fig. Here the label p denotes a generic label, such as p for NRQCD 
quarks, (n ■ p, p±) for collinear SCET quarks, etc. One residual momentum is removed for 
each sum over a label momentum. 7 Any remaining residual momenta appear as integrals 
in their own right, providing a proper implementation of the multipole expansion. 8 In 
using equations such as Eq. fllBD , one should imagine using Fig. || with explicit bins (i.e. 
hard cutoffs) for the k integrals. Once the final expressions are derived, one can evaluate 
the integrals in dimensional regularization, treating the bin sizes as infinite as discussed in 
section TT|. Eq. flip]) is only true if we sum over all p, including p = 0. The replacement 
in Eq. ( |IBD is what we will call the result for the naive integral, generically denoted with a 
tilde, I. 

We can finally formulate the zero-bin subtractions mentioned at the beginning of this 
article. The sum on p in Eq. ( |T6"D is over p ^ because the effective Lagrangian terms are 
a sum over p ^ 0. The restriction p ^ in Eq. (|l^) modifies the right hand side. In a 
Feynman graph, let F ({pi} , {ki}) be the integrand, including all the momentum conserving 
^-functions and label preserving Kronecker-5's at the vertices, where % runs over all the 
internal propagators. When integrating over a function F, the correct form of Eq. ( |l6f) is 
actually 

e / n ^ f > /n [ f - e p t ] • 

On the l.h.s. the sum is over all label momenta avoiding the zero bins, Pi = 0. On the r.h.s. 
we integrate pi over all of momentum space and the second term subtracts the contributions 
from regions j G U where one or more Pi vanish. The set of such regions U can be broken 
up into Ui where pi = 0, U# where Pi = and pj = 0, etc. The subtractions over U are 
defined iteratively by first subtracting over each Ui, then adding back U^, subtracting U^, 
etc. Note that since the Pi ^ constraint comes from the fields in the Lagrangian, it is 
implemented at the level of propagators in a graph, i.e. for each internal line, not for each 
loop momentum. Once the momentum conserving delta functions are accounted for the 
subtractions are implemented at the level of the full integrand. 9 



7 From the reparameterization invariance j58| in splitting p + k it is equivalent to think of this as first fixing 
the lattice of p's and adding the integrals over k, or as using the freedom in the choice of p's to fix k and 
then extending the sum over p's to an integral. 



8 See 1 59 for why this is relevant to power counting in NRQCD. 

9 If there are less sums over label momenta pi than integrals dki then these extra integrals over k's will 
appear on the RHS with corresponding dependence in F. These extra k's are momenta that are truly 
small for the physical process, see Ref. p4| . For simplicity this complication was suppressed in writing 
Eq. (|l7|), since it is not the most important aspect. 
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Each of the Pi ^ terms in Eq. (fP7|) represent the full label on a field, which for example 
will be a four- vector for soft gluons in NRQCD, g M , and the = q~n fJ- /2 + components 
for a collinear quark in SCETi. The proper subtraction integrand _p sub is obtained from F 
by assigning a scaling to all the pi appropriate for the zero-bin region, and expanding in 
powers of the momentum which vanishes in the zero-bin (as we scale towards this region 
with the power counting parameter in the sense of a standard OPE). For example, if we 
sum over g M 7^ for a soft momentum (i.e. order mv) in NRQCD, then we define F snh by 
an expansion of the integrand F in powers of q^ by assuming that q^ ~ v 2 , i.e. by assuming, 
for the purposes of the expansion, that q^ is usoft. This subtraction is done at the level 
of the integrand. The expansion is done to high enough order that the resulting integrand 
F — F suh vanishes in the zero-bin, i.e. vanishes as q^ — > 0. 

Note that there is a freedom to define a scheme which leaves a finite integrand in the 
scaling limit since this just moves finite pieces around between the matching and matrix ele- 
ments. We will use the scheme where all singularities are removed but not finite pieces. Also 
note that in a given subtraction some terms will be power divergences which in dimensional 
regularization are set to zero. This implies that if we had subtracted additional polynomial 
pieces they would not change the result of loop graphs in dimensional regularization since 
they integrate to zero. 

To show that the RHS of Eq. ( fPTD provides the proper implementation of the zero-bin 
independent of imposing hard momentum cutoffs we use a logic similar to section [ITT] in 
discussing Figs. Q-f|. The J2 p ^o can ^ e turned into a full integral if we add the p = bin, 
but we must subtract it again. The subtraction term has only an integral over residual 
momentum, but when we send the hard cutoffs on the sides of the zero-bin to 00 the 
subtraction is also integrated over all momentum. Thus we end up subtracting terms derived 
from the scaling limit of the original integrand integrated over all of momentum space. 
The full integrand on the r.h.s. of Eq. (\Tf\j ensures that we do not double count the zero- 
bin because the integrand vanishes when the loop momentum is sent towards the zero-bin 
momenta. This ensures there is no double counting in the IR. Any double counting in the 
UV is taken care of by Wilson coefficients and renormalization as usual. 

It is worth emphasizing that the result in Eq. (|T7|) applies equally well to the use of 
scaleless regulators like dimensional regularization, and to the case where hard Wilsonian 
cutoffs are applied to distinguish modes. For the Wilsonian case, consider the cutoffs as 
^-functions multiplying the integrand. In this situation the regulator ensures that the inte- 
grand is zero in the scaling limit so the subtraction terms all turn out to be zero, and the 
naive replacement in Eq. ( |T6| ) with the cutoffs gives the correct answer. 

The result in Eq. fliT| ) applies to NRQCD or SCET or any other quantum field theory of 
this type. It is necessary to avoid double counting the zero-bin momenta which correspond 
to different degrees of freedom in the effective theory. It provides a means for tiling the 
infrared regions of a quantum field theory with different degrees of freedom while avoiding 
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double counting. In the following sections we explore the difference between Eqs. ([16]) and 
(|T7|), and its consequences, with the help of several examples. 

Before leaving the general discussion it is worth emphasizing that Eq. ( |TTD together with 
the definition of the NRQCD and SCET degrees of freedom given in Figs. [5] and [| give 
complete coverage of all momentum regions where IR divergences can occur. We used 
multiple degrees of freedom to cover these IR regions because this facilitates setting up 
the proper EFT power counting expansion. In general one can look at combining regions 
together to describe a larger region with only a single degree of freedom. Doing so comes at 
the expense of making the power counting expansion difficult to formulate. In many cases it 
is actually unknown how to formulate the EFT expansion when regions are combined, thus 
necessitating multiple modes. We see that in general, the concepts of i) an EFT having a 
complete set of degrees of freedom to reproduce all IR divergences, and ii) the EFT having 
a valid power counting expansion, are tied together. We use this freedom to define NRQCD 
and SCET to cover the IR regions with our chosen degrees of freedom, so they reproduce the 
IR divergences. In this case proving that these EFT's are complete is equivalent to proving 
that their power counting expansions do not break down at any order. Demonstrating this is 
easier, since the power counting can only break down if we have missed a relevant operator 
at leading order in the expansion of some observable. All subleading operators are treated as 
insertions and do not upset the power counting. Thus, we see that constructing a complete 
EFT is equivalent to identifying the proper physical degrees of freedom in the leading order 
action, which is related to identifying the set of physical processes for which the EFT applies. 

V. ZERO-BIN SUBTRACTIONS IN NRQCD (NON-RELATIVISTIC PRO- 
CESSES) 

In this section we consider examples of the use of Eq. (|17|) for non-relativistic field theories. 
The results are quite general, applying whether the non-relativistic particles are quarks, 
nucleons, ions, or quasi-particles. The fields generating the potential are different in these 
cases, but the same momentum regions are important. To be definite we use a gauge theory, 
and so take our examples from non-relativistic QCD (NRQCD). 

A. NRQCD Soft Crossed-Box Graph 

It is helpful to consider a concrete example to study the consequences of Eq. ( |17D — we 
will start with the crossed-box graph in NRQCD. The full theory integral contains hard, 
soft, and usoft contributions, and we examine the soft crossed-box graph in the effective 
theory shown in Fig. |6]. Here pi and P2 are the external momenta, with r = P2 — pi the 
momentum transfer, and 2mE = = p|. The momentum of all the particles has been 
denoted as (label energy, residual energy; label momentum, residual momentum), and the 
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(0,E;Pi,0) (p ,E+k°;pi+P,k) (0,E;p 2 ,0) 



(p°,k ;p,k) (p°,k ;p-r,k) 



(0,E;-Pi,0) (p ,E + k°;p-p 2 ,k) (0,E;- P2 ,0) 

FIG. 6: Soft crossed-box graph in the effective theory. The zigzag lines are soft gluons, the double 
lines are soft quarks, and the single lines are potential quarks. For each line we show (label energy, 
residual energy; label momentum, residual momentum). 



external particles have been chosen to have zero residual momentum. 

In Feynman gauge the propagator for a soft gluon with momentum (p°, k°; p, k) is 
l/[(p ) 2 — p 2 + i0 + ], and for a soft quark with momentum (p°, k°; p, k) is l/[p° + ie]. We ne- 
glect overall factors such as color Casimirs and coupling constants, and focus on the integral 
for this graph: 

/cross V fdPk 1 1 1 1 

These soft propagators do not depend on the residual momentum components and so the 
propagator takes different values at each grid site in Fig. |3], but the same constant value for 
all the points in each box. Thus there is no change to the integrand when the label and 
residual momenta are combined into a continuous integration using Eq. which gives 

fcross f d °P 1 1 1 \ ( , q] 

s J (2Tr) D p° + iO+ p° + iO+ (p°) 2 -p 2 + i0+ (p ) 2 - (p - r) 2 + i0+ ' 1 ' 

However, Eq. (|T^) still includes the zero-bin contribution, which must be subtracted out. A 
more careful analysis taking account of the zero-bin for each soft particle propagator, and 
instead using Eq. dl7|) implies that the value of the soft crossed-box graph is not Eq. (|i9|) 
but rather 

rcross fcross rcross rcross l'On N \ 

l s — lg l x 1 2 [ZV) 

where the subtractions are: (a) 1% from the region (p° = 0, p = 0) and (b) I2 from the region 
(p° = 0, p = r): 



rcr 
1 1 



Oil f f 



{2tt) d p° + i0+ p° + i0+ {p ) 2 - p 2 + i0+ -(r) 2 + ie ' 



jc^oss = / _d^^ 1 1 1 (21) 

2 ' (2tt) d p° + i0+ p° + i0+ -r 2 + z0+ (p ) 2 - (p - r) 2 + z0+ ' K J 
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(a) (b) 

FIG. 7: Ultrasoft corrections to potential scattering related to the crossed box graph. 



By the shift symmetry in p we have /^ ross = /£ ross . 

The /cross su btraction comes from the region where the (p°, k°; p, k) gluon is usoft. Sim- 
ilarly, the /| ross subtraction comes from the region where the other gluon becomes usoft. 
Subtractions from the region where the quark is potential (p° = 0) vanish, as do the double- 
subtractions where the regions for overlap with p° = 0. This is because for the crossed 
box all the p° poles are on the same side of the contour of integration. 

The subtraction i^ ross avoids double counting the usoft graph shown in Fig. |7|(a) and 
similarly /£ ross avoids double counting the usoft graph Fig. |7|(b). The usoft graphs depend 
on external usoft variables such as the energy. If these are set to zero, then the integral 
Fig. |7|(a) is J^ ross ; including the omitted color factors. The reason that /™ oss does not depend 
on external usoft variables, but Fig. [F](a) could, is because /™ oss is obtained by considering a 
soft graph, and then taking its usoft limit. The effective field theory Feynman rules require 
that all usoft momentum be expanded out while considering soft diagrams. Thus I^ oss is 
the same as Fig. 0(a) with the external usoft variables expanded out. 

Prior to the subtraction the soft crossed-box integral in dimensional regularization is 



4 7r 2 r 2 



— + ln (*T 



(22) 



The total subtraction /= ross + /cross gj ves 



rcross i rcross 



4 7r 2 r 2 



(23) 



so the final result for the crossed-box integral is 



rcross rcross jc 
1 S ~ i l ~ 1 2 



An 2 r 2 



In 



(24) 



The subtractions have converted the 1/erR divergence in J^ TOSS into a 1/euv ultraviolet di- 



vergence. In Ref. 0, 0, it was argued that 1/em divergences in soft graphs should be 
converted to ultraviolet divergences by a pullup mechanism and included in the computation 
of anomalous dimensions. We see that the zero-bin subtraction automatically implements 
this conversion. An important feature for NRQCD is that the soft graph defined with the 
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(0,E;p,,0) (p ,E+k ;p 1+ p,k) (0,E;p 2 ,0) 



(pVlP.k) g 



1 (P°,k°;p-r,k) 

3 



(0,E;-P„0) (-p°,E-kVp,-p.-k) (0,E;-p»0) 



FIG. 8: Soft box graph in the effective theory. The zigzag lines are soft gluons, the double lines 
are soft quarks, and the single lines are potential quarks. For each line we show (label energy, 
residual energy; label momentum, residual momentum). 



zero-bin subtraction is infrared finite, and has a well-defined renormalized value independent 
of any cancellation with usoft graphs. 

A similar conversion from infrared to ultraviolet divergences also occurs for collinear 
graphs in SCETt as we show in section |VI[ 



B. NRQCD Box Graph 

The zero-bin subtraction has another important consequence — it gets rid of pinch sin- 
gularities. Consider the soft box graph in NRQCD, shown in Fig. |8|, and follow the same 
procedure as for the crossed-box graph. The only difference from the crossed-box is the 
replacement 

1111 

(25) 



p° + iO+ p° + iO+ p° + iO+ -p° + iO+ 
in Eqs. (JLl^pi~|) due to the change in momentum routing through the antiquark line. The 



integrals analogous to those in Eqs. (|i9|j2l|) are 

d D p 1 1 1 



I 



box 



s ~" J (2n) D p° + m -p° + m (p ) 2 - p 2 + i0+ (p ) 2 - (p - r) 2 + zCF 
d D p 1 1 1 1 



rbox 



(27r) D p° + iO+ -p° + iO+ (p°) 2 -p 2 + i0+ -{r) 2 + ie 



/box = f d D P 1 1 _A 1 (26) 

2 J (2n) D p + iO+ -p° + iO+ -r 2 + i0+ (p ) 2 - (p - r) 2 + i0+ ' 1 ; 

where the usoft subtractions l\°2 are for (p° = 0, p = 0), (p° = 0, p = r) respectively. In 
addition one also has a potential subtraction for p° = 0, 

/b ox = [Pv__!: 1 1 1 m 

3 J (2rr) D p° + iO+ -p° + i0+ -p 2 + z0+ -(p - r) 2 + z0+ ' V ; 
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Now this p° = subtraction overlaps with the usoft subtractions, so we have to add back 
the double subtractions, the (p° = 0) limit of Jj 5 ™: 



rbox 



rbox 



1 



1 



1 



(2ir) D p° + i0+ -p° + i0+ -p 2 + i0+ -r 2 + i0~ 
d D p 1 1 1 1 



(2tt) d p° + i0+ -p° + i0+ -r 2 + i0+ 
The complete expression for the soft box graph is 



IP 



+ i0~ 



(28) 



rbox 
J 5 



rbox 
J 5 



box 



rbox 
J 2 



rbox _j_ rbox _i_ rbox 



(29) 



Both 7| ox 



and /I"™ 



have pinch singularities in the p° integral, from the poles at p° = ±i0 + , 
and are ill-defined. However, for the result in the effective theory, we don't need the separate 
integrals, but only the difference /| ox — I^ ox , which has no pinch. We have 



rbox 



rbox 



d D 



P 



(2ir) D p° 



i0~ 



X 



(p0)2 _ p 2 + t Q+ (p )2 - (p - r ) 



i0- 



i0+ (p 



^0^ 



30) 



One can evaluate the p° integral in Eq. Q3Tj) using contour integration. The result is the same 
as doing Eq. ( [Of) by contours and dropping the pinch pole at p° = 0, since the integrand of 
Eq. (|30D has no p° pole and the subtraction term does not introduce new poles in p°. This 
prescription for the soft box graph is what was used in Refs. || [Ll], 0, [49], ^0f, but we 
now see how the effective theory automatically gives this result. 10 

The double subtractions and Jg™ remove the pinch poles at p° = for the sub- 
tractions Ii° x and J2° x respectively, so that J] 3 ™ — / 4 box and /j™ — -^5° x are f ree of pinch 
singularities. This justifies ignoring the p° = pole in the calculation of these integrals. 
The l\°2 usoft subtractions convert the infrared divergences in /| ox into ultraviolet diver- 
gences, just as they did for the crossed-box. Prior to the subtraction the soft box integral 
in dimensional regularization is 



rbox 



rbox 



47r 2 r 2 



The total subtraction is 

rbox 



IT 



rbox I rbox 
J 4 1 2 



rbox 



1 , (V? 
— + ln ^ 

em V r 2 



1 

em 



4 7r 2 r 2 



euv 



(31) 



(32) 



At one-loop there are many prescriptions that lead to the same result, that the pinch is dropped in the con- 
tour integration. Examples include split dimensional regularization and the principal value prescription. 
The formula in Eq. (0) can be applied at any order, and makes adopting a prescription moot. 
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so the final result for the box integral is 



rbox rbox rbox rbox rbox , rbox , rbox 

1 S — 1 S ~ i l ~ 1 2 ~ J 3 ' A "r J ; 



euv V r 2 



(33) 



The subtractions have converted the 1/eiR divergence in /| ox into a 1/euv ultraviolet diver- 
gence just like the cross-box. The zero-bin subtraction has removed the infrared divergences 
and the pinch singularities, since these regions are properly taken care of by usoft and po- 
tential graphs, respectively. The properly defined soft box graph is infrared finite, and has 
no pinch singularity. 

The use of Eq. ( |TT| ) with the zero-bin subtraction works at higher orders as well. One 
can check explicitly that the subtractions remove the pinch singularities in the double box. 
If the loop momenta for the two single-box subgraphs are called p and £, the subtracted 
double-box is given by subtracting the region where p° = and where £° = 0, and adding 
back the region where p° = £° = 0. This gives an expression for the double box which is free 
of both single and double pinch singularities. We have checked that the subtracted three 
gluon exchange graphs give the correct contribution to the two-loop static potential |H| in 
Feynman gauge. 

The standard computation of the Coulomb potential at 0(a 2 ) in QCD is free of IR 
singularities for a different reason, because the IR divergences in the box and crossed-box 
cancel against the vertex and wavefunction diagrams. With the zero-bin subtractions, the 
scaleless soft vertex and wavefunction diagrams are set to zero, and the box and cross-box 
together with the non-Abelian vacuum polarization and Y-graphs give the complete IR finite 
answer. At this order there is also no overall usoft contribution to this four point function, 
and the same is true for the Coulomb potential at two loop order |JT], |62[. At three loops 
the zero-bin subtraction removes the ADM singularity in the Coulomb potential |15|, |64| . 



A one-loop example at 0(v 2 ) where the subtractions do not cancel in the sum of diagrams 
is discussed in the next section. 

C. Results for the Box and Crossed-Box at order 0{v 2 ) 

In this section, we study the soft box and crossed-box graphs to second order in the v 
expansion, i.e. to the same order as the spin-orbit, Darwin and tensor-force contributions to 
the QQ potential. At order v 2 , the naive soft box and crossed-box have IR divergences and 
there are also non-trivial contributions from usoft diagrams. We will summarize results for 
these graphs to illustrate how the zero-bin subtractions work. This example also illustrates 
a case where Fj ub in Eq. flTT| ) involves a series of terms. 

Prior to any subtractions, the necessary diagrams are simply given by all quark- ant iquark 
scattering diagrams that are derived using the HQET Lagrangian up to 1/m 2 . The full 



integrands are lengthy and we refer the reader to Refs. R 42]. After using standard tricks 
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to trivialize the numerator momenta we are left with the basic integrals 

d D k 1 1 



J(a,(3) 



E 

p<=Z c 



(2n) D ( P +iO+) a (- P °+iO+f [Go°) 2 -p 2 +i0+] [(p0)2_( p _ r )2 +i0 +]' 

(34) 

where Z c = {p° ^ 0,p^ ^ 0,p M ^ (0, r)}. The subtractions that account for p° ^ and 
remove the pinch singularity from the first two denominators involve 5 derivatives of the 
second two denominators where 5 is the nearest integer < (a + (3). For this particular 
computation this is equivalent to ignoring these poles in the contour integral. 

As in the previous section the removal of the remaining constraints is similar, with or 
without the pinches, so we will consider the case (5 = for simplicity. The naive integral 
and its subtractions for {p^ ^ (0, r),p^ ^ 0} respectively are 

d D p 1 



J(a,0) 
Ji(a,0) 

M<x,0) 



(2ir) D (p° 
d D p 



io+) a [(p°) 2 -p 2 + io+] [(p 1 



0\2 



(2tt) d (p° + io+y 

d D p 



,0)2_( p _ r )2 + i0 H 



p 

a-2 r 

E 1 

fc=0 



-r) 2 + i0+] ' 
-2 (p - r) ■ r] 



1 y [-2 P -r] fc 

(2n) D (p + i0+) a [(p°) 2 -p 2 + i0+] f^ Q [-r 2 ]^ 1 ' 



(35) 



where in J\ and J2 we have dropped terms that are obviously zero. Note that here removing 
the zero-bin requires a series of subtractions obtained from expanding the naive integrand 
about the zero-bin values. By translation invariance in dimensional regularization the two 
subtraction integrals are equal, Ji(a,0) = J2(a,0). 

For a — 1, the subtractions are zero, J\ = J2 = 0, and J(l, 0) is finite. For any other odd 
a, both the naive integral and subtractions give zero. For even a < 0, J(a, 0) is UV divergent 
and the subtractions give zero. Finally for even a > 2 the base integral is IR divergent and 
the subtractions convert this to a UV divergence for J(a, 0) = J(a, 0) — J\(a, 0) — J2{&, 0). 
As an example, consider a — 4. The naive integral is 

1 



J(4 - 0) = i^ 

and the subtractions J\ = J2 are given by 

d D p 1 



16 \ 
3rV 



+ ln 



(36) 



^(4,0) 



(2tt) d (p° + i0+) 4 [{p°) 2 -p 2 + i0+] 
i ( 16 \ r 1 1 



1 2pr 4(p-r) s 



= + + 



16tt 2 



3r 4 7 



Thus the full integral is J(4, 0) = J(4, 0) - 2 J 2 (4, 0) giving 



J(4,0) 



16tt 2 



16 \ 
3rV 



+ In 



+ 2 



(37) 



(38) 
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which does not have an IR pole. The same is true for all even a > 2. 

Lets consider the sum of all order v 2 NRQCD soft exchange diagrams (boxes, cross-boxes, 
and triple gluon graphs), and the vertex and wavefunction graphs, all computed in Feynman 
gauge with equal mass quarks and antiquarks. Using the naive J integrands we find 11 



5 _ 



Ci(l®l)-^(T A (8)T A ) 



1 4 

+ 



+ C A {T A ®T A ) 



4e 



13 1 

+ 



uv 



+ A 



1 



2euv e IR 



3e m 
11 



(P 2 + P /2 
2r 2 



(39) 



14 \ 



18e uv 3e 



S. 



vertex+w.fn. 



1 -^{T A ®T A ) 



3V 2 



1G, 



IR 



P 2 + P /2 

2r 2 



3euv 3eiR 
1 1 



)) 



36e uv 6e m 

S 2 3A T 
"3 2 ~ 12 



x 



1 



1 



CUV e IR 



for the pole structure. We used the notation in Ref. |§ where in SU(N C ) the color coefficients 
d = (iV 2 - 1)/(4N 2 ), C d = N c - 4/JV c , C A = JV C , C F = (iV 2 - 1)/(2AQ, there are two color 
structures (1 <S> 1) and (T A eg) T A ), and A, S 2 , and T are spin and momentum dependent 
structures, 



(J\ + a 2 



A = -i 



,S-(p'xp) 



3 r ■ <7i r ■ cr 2 



(40) 



The sum of diagrams with naive integrands is 



S = l ^±{ Cl {l^ 1! 

77T 



4 \ 



+ C A (T A ®T A 
7 
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)+(T 
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(J-J-)-o,( J- 
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1 \ , (p 2 + p' 



2r 2 



3e uv 3e 



IR 



+ A (2e uv ) + S (l8e uv ) + T (36e uv ). } 



3e 
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(41) 



After subtracting the zero-bin contributions to get the proper integrals J(a,/3) we find 



7 



3euv 



C d (T A ® T A ) 



7 



12e uv 



C A (T A ® T 
5 



f— ) 

V12e uv / 



2r 2 V3euv^ \2euv' V18e uv / V36e uv /J J v ; 



For the purpose of this example we have set Wilson coefficients in the HQET Lagrangian to their tree 
level values. For the summation of logs in Refs. (4^, ^ J53|, their renormalization group evolution of 
course had to be kept. 
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At 0(v 2 ) there are also UV divergences from the usoft graphs which can be found from 
Ref. [j3S| . Using dimensional regularization for both the UV and the IR they give: 

, = 4(i.I)M( 181) + [ Cj 4A«!^]p- gJ -)}. 
m z Veuv cir/ I 3 L 3 3 3 2r z J J 

(43) 

We see explicitly that the usoft graphs have UV divergences which match up with the fake 
IR divergences from the unsubtracted soft computation in S. The true soft computation 
gives an IR finite result and the usoft contribution exactly matches the IR divergences in 
the full theory computation, see Ref. ||. To interpret the UV divergence in Eq. ([|3|) as 
occurring at the hard scale, it is crucial to make the zero-bin subtractions to avoid double 
counting in the soft region. 

VI. ZERO-BIN SUBTRACTIONS IN SCETi (INCLUSIVE PROCESSES) 

SCET is another theory with correlated scales, and with multiple fields for the same 
particle. As a result, one expects the zero-bin subtraction to also apply in this theory 
(yielding a pullup here too). In this part, we consider examples of the application of Eq. ( |TTD 
to SCETi which has collinear fields appropriate for the description of perturbative energetic 
jets, and non-perturbative usoft fields. 12 The new feature of SCETi is the appearance of 
double-logarithmic divergences at one-loop order. In the following discussion we show how 
the zero-bin subtraction works in this case. 

We first consider the heavy to light vertex diagram in Fig. ^ which appears in processes 
such as the b — > 57 transition magnetic moment operator needed to compute inclusive 



B — > X s 7 decays. In section [VI A| we work on-shell with finite cutoffs and demonstrate 
that the infrared divergences of the full theory are reproduced in SCETi. The subtraction 
from Eq. (|P7D in the collinear diagram is required for this to be true. To demonstrate that 



Eq. (|T7|) is independent of the choice of ultraviolet and infrared regulator, in section [VI B 



we 



consider the more standard choice of dimensional regularization with an offshellness infrared 
regulator, and explain how the zero-bin subtraction works for this case. In section |V1 C 



we 



treat the example of the current relevant for inclusive two-jet production where we have 
collinear fields in two directions, n and n, and usoft fields. The Yl P7 to are a ^ so important at 
tree level and for phase space integration as demonstrated by the 7* — > qqg example that 
we take up in section |V1D . 



These collinear modes are sometimes called hard-collinear, and the usoft modes are then called soft. 
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FIG. 9: Heavy-to-light vertex diagrams, a) full theory, b) SCETi collinear graph, and c) SCETi 
usoft graph. 



A. Onshell Integrals and a Cutoff Regulator for B — > X s j 

For b — > S7 at lowest order in Aqcd/-^ we have the SCET current || 

J® =C(u)[(£ n W) u Th v ]. (44) 
After making the decoupling field redefinition Q on the collinear fields this becomes 

=CH[(^).r(rt^)]. (45) 



We start by making use of the current in Eq. ( f44"|) and will discuss the equivalence of using 
(|45|) at the end of this section. The incoming heavy quark momentum is p^ = mv^ with 
v 2 = 1 and the outgoing light quark momentum is pf" = p^n^ /2 where n 2 = and p- = n-p. 
Both the incoming and outgoing quarks are taken onshell. 

In this section we first demonstrate the effect of the pullup on collinear diagrams, prior 
to making a specific choice of regulator. We then use infrared cutoffs Q± and f2_ on p± and 
p_ in both the full and effective theories, so that the loop momenta are restricted to the 
region q\ > Q 2 ± and q 2 _ > Sl 2 _. For the usoft graphs, we use a ultraviolet cutoff A_, while for 
the collinear graphs, we use an ultraviolet cutoff Aj_. Hard cutoffs make the computation 
of anomalous dimensions more difficult and in more generic diagrams would require gauge 
violating (and power counting violating) counterterms to restore these symmetries. Our 
focus is on showing how the IR divergences are reproduced for a particular example where 
these problems do not occur, so for the purpose of this computation these issues are not a 
concern. 

The part of the full theory diagram in Fig. ||a with the double logarithmic infrared 
divergence involves the integral 

i*-"7 = \1l± ^1 (46) 

ful1 J (2-K) D (q 2 + m+)(q 2 + 2p b -q + M+)(q 2 + 2p-q + i$+)' 1 ' 

which is ultraviolet finite. Taking pj, = rrib(n + n)/2 and < p~ < mf,, we use the identity 
d D q = d(n ■ q)d{n ■ q)d n q±/2 to write the measure in light-cone variables, where here the 
exponent n = D — 2 is the dimension of the J_-space, not to be confused with the light-like 
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vector n^. Next we perform the n-q integral by contours. There are poles for 

9± - n + • \ (gi+m 6 n-g) . . 

n-g = — «U sign(n-g) , n-g = = i\j signm- q+rrib) , 

n-g n-q+rrib 

q 2 

n-q——- i0 + sign( n-q+n-p) , (47) 

n-q+n-p 

which gives three poles above the axis for q~ < —nib, one moving below in the region 
— nib < Q~ < —p~, two below and one above for — p~ < q~ < 0, and all below for q~ > 0. 
Thus only the middle two regions contribute. We will drop the integral over the interval 
— m < n-q < — p_, since it is finite in the UV and IR. This gives 

-ff.^T 7 = — / dn-q - — - — t — ~ ; , on . . + finite , (48) 

fuU 2nJ_ p _ (27r)"[ql + (n-qy\(ql) 1 } 

The usoft graph in SCETi is shown in Fig. |9]c and just has an integral over residual 
momentum. It is important to recall that the leading order SCET Lagrangian involves 
a momentum space multipole expansion for the residual momentum 0, so that only the 
residual n ■ k momentum appears in the collinear quark propagator. The integral is 

it-* = rsi^ ^ = i av* dnk ' n - k 



us J (2n) D (k 2 + iO+)(vk + iO+)(n-k + iO+) 2tt J-ao (2tt)" [k 2 ± + (n-k) 2 }(kl ) ' 

(49) 

For the collinear graph, Fig. |||b, we have the label loop momentum (p = n-qn^/2 + q^_ and 
residual loop momentum g£\ and we will denote q^ = g M + n-q r n^ /2. The original integral is 

/6 _ S7= >p f d D q r 2n-{q + p) 

° mtt-p J ( 27 r) D (n-q + i0+)(q 2 + 2p-q + i0+)(q 2 + i0+) ' 1 ' 

Eq. ( pTT| ) is used to take into account the subtractions from the zero-bins. For q ^ we 
examine the scaling (p ~ A 2 , for which case the loop- measure scales as [A 8 ] and combines with 
the integrand to give: [A 8 ]/[(A 2 )(A 2 )(A 4 )] ~ A , so there is a non-trivial subtraction for this 
region. For q ^ —p we examine q~ + p~ ~ A 2 , g -1 ~ A 2 and have: [A 8 ]/[(A°)(A 4 )(A 2 )] ~ A 2 , 
so this zero-bin can be ignored. When we combine the sum over label momentum with 
the integral over residual momentum we get the naive result Iq and a subtraction Jo from 
Eq. ©: 

fb->ry f dD( l 2n-(q + p) i f° d n q ± n-(q + p) 



(2n) D (n-q + i0+)(q 2 + 2p-q + i0+)(q 2 + i0+) 2tt J_ p - * (2ir) n (n-q)(n-p)(q 2 ± ) ' 

f~<n = [_f± 2 ±l = ± f° dn a n 'P f 51 ) 

J (2n) D (n-q + i0+)(n-qn-p + i0+)(q 2 + i0+) 2ix J _^ H (2tt)™ (n-q)(n-p)(q 2 ± ) ' V ; 

Here the zero-bin subtraction I^ S1 is obtained from the — > QX 2 scaling limit of the J^p S7 
integrand and avoids double counting for the usoft region of momentum space. To double 
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logarithmic accuracy the n ■ q in the numerator of f^ sl can be dropped so 
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(52) 



We see that the subtraction integral changes an infrared divergence in Iq at n-q = into a 
ultraviolet divergence for n-q —>■ — oo. Since we can see this at the level of the integrand it 
is obviously independent of the choice of ultraviolet and infrared regulators. 

With the prescribed cutoff regulators and _L-spacetime dimension n = 2, these integrals 
can be evaluated to give 
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(53) 



The full result for the collinear graph is therefore 
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(54) 



and we see that the zero-bin subtraction l h ^ sl has converted an IR divergence ln(0_) for 
the q~ variable in I^ S1 into a UV divergence, ln(A_). The sum of the SCETi effective 
theory contributions gives 



6— >S7 
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+ ln 2 f^)-ln 2 (^ 
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(55) 



The first two terms on the r.h.s. contain the infrared divergences and exactly reproduce 
these divergences in the full theory result Z^ 57 . Furthermore, the last two terms in Eq. ([55]) 
depend only on the ultraviolet cutoffs and the large label momentum p~ and can be compen- 
sated by a counterterm for the current in SCETj. If Jq^ S7 in Eq. (|54|) had been left out, then 
we would not properly reproduce the IR divergences in the full theory result. Furthermore, 
without Jq^ S7 , the ultraviolet cutoff dependent term would have cross terms ln(A_) ln(f2j_) 
and ln(Aj_) ln(fi_) and it would not be possible to cancel the cutoff dependence by a coun- 
terterm independent of the IR regulator. 

The above calculation was performed for the current in Eq. fl44]) . Since our regulator 
leaves all external lines onshell we obtain exactly the same results if we had started with 
the current J'^ in Eq. (fi5|) , which is obtained after making a field redefinition involving the 
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Wilson line Y. Since we work onshell the two forms of the current are equivalent, and the 
Feynman rule from the Wilson line Y give exactly the same integral in Eq. (|49|). Thus our 
implementation of a cutoff IR regulator does not destroy the eikonal factorization embodied 
by the field redefinitions involving the Wilson line Y. This property of the field theory is 
not maintained with the offshellness IR regulator which we consider in the next section. 
This should be considered as a fault of this IR regulator as pointed out in Ref. p5] . In 
Ref. [BS] an energy dependent gluon mass regulator was studied which also preserves the 
field redefinition. 13 



B. Offshell Regulator with Dimensional Regularization for B —> X s j 

We now repeat the calculation of the effective theory diagrams in the previous section but 
keep p 2 7^ to regulate the infrared and use dimensional regularization for the ultraviolet, 
D = 4 — 2e. The full theory integral is 

m " J (2n) D (qi + i0+)(q2 + 2p b -q + i0+)l(q + p)'> + i0+] ' 1 ' 

The SCET integrals are 

jb—^s'y 



d u k 1 
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(2tt) d (k 2 + iO+)(v-k + iO+){n-k + p 2 /n-p + i0+) 
d D q 2n-(q + p) 

(2tt) d (n-q + iO + )[(q + p) 2 + i0 + ](q 2 + i0+) ' 



J (2vr) D (n-q + i0 + )(n-qn-p + p 2 + i0 + )(q 2 + i0 + ) ' 1 ' 

Again, one can see that as n ■ q — > the difference 7^ S7 — Jq~^ S7 does not have an infrared 
divergence from this region. However in Ic alone, there is an infrared divergence from this 
region that is not regulated by p 2 ^ 0. It is regulated by dimensional regularization, and so 



13 Ref. |5f| also argued that the n-k — > divergence must be treated as a UV in the EFT since it comes from 
angles opposite to the collinear direction. The renormalizability properties of field theory only appear for 
large momenta, and the zero-bin turns this divergence into a true UV divergence. One must be careful 
about the distinction between angles for particle and antiparticle poles when determining that the n-k — * 
divergence is IR. 
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contributes to the 1/e singular terms. Evaluating the above integrals we find, 
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where we have distinguished between ultraviolet and infrared divergences. Here we see that 
the zero-bin contribution Iq~^ sj is responsible for canceling IR divergences in I^ SJ that were 



not regulated by the offshellness, 
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Therefore the sum of the SCETi contributions gives 
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The last term reproduces the infrared structure of the full theory result. The first two terms 
are canceled by a counterterm. The third term contributes a finite contribution to the hard 
Wilson coefficient of the heavy-to-light current in matching onto the full theory. Again we 
see that the contribution from I^ S1 is necessary in order for the infrared divergences in the 
full and effective theories to match up. The ellipses in Eq. ( |6~0"1) denote 1/e, single log, and 
finite terms that we have not bothered to display in the quoted results, but which have the 
same desired properties. 

The computation of SCET anomalous dimensions in Ref. @ and all subsequent papers 
used the entire 1/e divergent terms in J^ S7 to compute the anomalous dimension. As we 
have shown above, some of these divergences are, in fact, infrared divergences. The pullup 
mechanism is needed to convert these into ultraviolet divergences which can then be canceled 
by local counterterms in the effective theory, and so properly contribute to the anomalous 
dimension. 

One can see the problem with having J^p S7 and no subtraction term in another way. 
Consider adding a very small gluon mass to the collinear calculation, where m 2 <C p 2 . 
The gluon mass can only affect the IR, and we should obtain the same form for the 1/euv 
divergences if we expand in m 2 before or after the integration. If we consider expanding after 
the integration, then performing the q + integral by contours followed by the q± integration 



34 



gives 
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Thus with m 2 ^ the IR singularity at q 



is no longer regulated by dimensional 



regularization in I^7* S1 or Iq 



however the difference I, 
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remains well 



defined. Here the 1/e terms in I C ~* S7 do not give the correct counterterm structure even if 
we set €ir = euv- Letting q~ = —xp~ gives 
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and in the limit m 2 <C p 2 this reproduces Eq. fl5T?|). Thus with the zerobin subtractions the 
1/euv divergences in the effective theory are independent of the choice of IR regulator. 

Note that in this section it was crucial to use the current in order that taking 
P 2 7^ provides the same IR regulator in the full and effective theories. For the SCET 
current J'^ this is no longer possible. Working with this current, only onshell IR regulators 
should be considered. This happens because the field redefinitions involving Y 7 s modify the 
LO collinear Lagrangian, rather than just subleading terms, and are therefore sensitive to 
regulation of the propagator. 



C. Production of n-n jets 

As another example of the zero-bin subtractions in SCETi we consider the one-loop 
diagrams contributing to two jet production, 7* — > qq. The degrees of freedom required 
in SCETi are those pictured in Fig. |]a. In this case we have subtractions for both the n- 
collinear and n-collinear fields, which ensure that they do not overlap with the usoft region. 
The leading order SCETi current is |6^ 

= c{u,u')^ n w n ) u r{wknh . (63) 

If we make the decoupling field redefinition || which encodes the eikonal coupling to all 
collinear quarks and gluons then becomes J'^ = C(u,u')(^ n W n )^Y^Y n 'y^{W^n) LlJ '. 
Below we work with the current since we will use an offshellness IR regulator. 

The one-loop vertex graphs are shown in Fig. pi]. Wavefunction graphs are not shown, but 
in Feynman gauge the collinear gluon wavefunction renormalization for a collinear quark is 
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FIG. 10: Vertex diagrams for n-n production, a) full theory, b) SCETi collinear graphs, and c) 
SCETi usoft graph. 



equal to the full theory result and the usoft gluon contribution vanishes. Working to double 
logarithmic order and breaking up the collinear terms into the naive result and subtractions 
the relevant integrals are 
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(64) 



We have kept offshellnesses, p 2 7^ and p 2 7^ to regulate the IR and will use dimensional 
regularization for the UV. Much like the heavy-to-light computation, this does not regulate 
all the IR divergences in the naive collinear integrands, I^ rod and /| >rod . Note that we took 
n-{p + q) — > n-p in the numerator of the collinear graphs since we only examine the double 
logarithms and have made a corresponding approximation in If^ d - Evaluating the full 
theory integral we find 14 
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(65) 



The usoft loop graph in the effective theory gives 
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14 In the computation of the cross-section in the full theory the ln(p 2 ) IR divergences are canceled by 
analogous IR divergences in the bremsstrahlung graphs. 
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For the n-collinear naive integral and subtraction we find 



jprod 



rprod 
'raO 



167T 2 

i 

16tt 2 



eiRfuv ClR 
2 2 \ r 1 



euv em ' l> euv 



— p' 
+ In 



-In 2 



-p^ 



— p* 
-In 



+ 

n-p 



em euv 



In 



n-p 



+ 



so the full n-collinear result is 
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(68) 



Just as in the b — > 57 example the subtraction terms remove the l/er R poles, and the IR 
in the complete collinear integral is regulated by the offshellness. The ellipses denote 1/e, 
single log, and finite terms that we have not bothered to display in the quoted results. The 
results for the n-collinear terms are similar 
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Adding up the SCETt integrals, J s p c r e ° t d = I^°J t + J? rod + JP rod , we find 
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(70) 



In the last line, the first two terms are removed by a counterterm in MS and the third term 
contributes to the Wilson coefficient C prod in the one-loop matching, see Ref. JOJ. The 
fourth term exactly reproduces the IR divergences in the full theory result. 

In the above computation there was an interplay between the usoft loop and the n and 
n collinear loops which combine to reproduce the IR of the full theory. The exact way in 
which these IR divergences combine depends on the choice of IR regulator as we saw in the 
b —>■ 57 example. Again we see that the zero-bin subtractions are important to correctly 
reproduce the IR divergences once we distinguish between euv and ei R . 

D. 7* — > qqg : The Zero-bin at Tree Level and for Phase Space Integrals 



In this section we show how the zero-bin is kept track of in tree level computations. 
One has to consider it even at tree level because in SCET multiple fields are present for 
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FIG. 11: Regions R=A-G for 7* —* qq g where the particles become usoft and collinear. For regions 
A,B,C the graphs in SCET are shown. The soft regions are slightly exaggerated for visibility. 



the same physical particle. We will demonstrate how double counting is avoided in fully 
differential cross sections and how the zero-bin subtraction affects phase space integrations. 
In the context of proving factorization in Drell-Yan, subtractions which avoid overcounting 
in phase space regions have been considered in Ref. [p7| . 

Consider the high energy process j*(q) — > q(pi)q(p2)g(P3) in the rest frame of the 7* 
with q 2 = Q 2 ^> Aq CD . This is a basic ingredient in two-jet and three-jet production, which 
were considered in SCET in Refs. [7[J. We take the full theory production current, 



J = ip^^ip. The external lines have p 2 = and we define dimensionless momentum fractions 
Xi = 2q ■ Pi/q 2 so that momentum conservation reads 2 = x\ + x 2 + £3- Computing the 
phase space integrals in (4 — 2e)-dimensions, the standard full theory cross section from the 
bremsstrahlung graphs is 



1 rfcr full CpOi s fi 2t 1 x\ + x\ — e(2 — x\ — x<i" 2 



cr dx l dx 2 2tt q 2e r(l-e)(a;i + X2-l) e (l-a?i) 1+e (l-x 2 ] 



l+e 



(71) 



where do is the Born cross section in dimensional regularization, a = (47ra 2 /Q 2 ) Y2f e< f + 
0(e) with a sum over the quark charges e/ = 2/3 or —1/3. 

Taking x\ and x 2 as the independent variables we have the phase space shown in Fig. |TT|. 
In SCET the different regions of the phase space plot are described by distinct EFT diagrams 
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with the particles being created by either a collinear or usoft field. The relevant limits are 
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(72) 



where 2 denotes generic values not near the two ends. In Fig. |TT]the SCET graphs for regions 
A,B,C are shown. We will compute the sum of the square of SCET diagrams for each region 
and compare them with the full theory result for the double differential cross section and 
a single differential cross section. In the 7* rest frame with 2q ti = Qn^ + Qn^ there is still 
a rotational freedom in the _L-plane which we can fix in performing the calculations. For 
A,B,C we take p^ = 0, while for E,F it is more convenient to take = 0. 

Computing these SCET graphs with the phase space integrals in dimensional regulariza- 
tion we find cross sections in each of the regions R 
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Here Ar is the amplitude in region R divided by Z 
diagrams we find 
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and 



, 1 + (X! - X 2 f 
\ A D\ X i+ X ^2~ 7T— + <-A e J 



X1+X2 7^ 2 

xi,a + 



1 2 



|2 



1-e 



SSo^ 2 (1-^)' 



l + x\ e (1— Xi) 



G\ 



|2 



11 / 1— X2 
XI + 1 

x 2 y^O 



xi+x 2 7^ 1 
^1,2 7^ 1 



(l-xi)(l-x 2 ) (l-x 2 ) 

0(e) . (75) 



(1-Xi)(l-X 2 ) 



Here « n , « s , and u are an n-collinear spinor, an n-collinear spinor, and an usoft spinor re- 
spectively, all with relativistic normalization. For the regions A, B, C, we explicitly show the 



amplitudes that follow from the SCET Feynman diagrams in Fig. [LI] and include e-dependent 
terms in the results. The amplitudes for A and B follow from the LO Lagrangians and LO 
SCET production current in Eq. (|63|). The result for C requires an insertion of the sublead- 



ing Lagrangian C^j = q us W^ig ^^n+^-C- J7T], where the field strength igTfi]_ = [in-D^, ip±\- 
In Eq. (p^) we have translated the zero-bin restrictions on the large momenta of collinear 
particles to restrictions on x\ and x 2 as shown on the RHS of the equations. 

It should be obvious from the form of Aa~Ag that one can not simply add the SCET 
diagrams to reproduce the doubly differential cross section in Eq. (f7l|) . The point is that 
the effective theory results do not overlap, as made explicit by the sums which exclude the 
zero-bins, and constrain the valid region of phase space. Given an x\ and x 2 , only one of the 
effective theory expressions is relevant. It is straightforward to determine which one once 
we specify parametric definitions of the scaling limits in Eq. (|72|), and pick values for X\ and 
x 2 . 15 The SCET diagrams in this region reproduce the full theory double differential cross 
section order by order in the expansion. Thus, it is crucial to take the zero-bin into account 
even at tree level in order to avoid double counting. 

Often we would like to deal with a less differential cross section which involves integrating 
over kinematic variables. In this case we should implement the zero-bin subtractions in the 
phase space integrals using Eq. ([H]) to avoid double counting when combining regions. As an 
example of the zero-bin subtractions in phase space integrals we consider the 7* — ► qqg cross 
section da jdx\ for fixed x± — 1 — S with S ~ A 2 and nonzero. In the full theory the single 
differential cross section for x\ — > 1 is obtained by integrating Eq. ( fTTJ ) over 1 — x\ < x 2 < 1 



15 The particles in the final state can be treated as observed, by a measurement of the final state. For 
each final state particle, one can assign a label p and residual momentum k as given by the binning of 
momentum space, and classify particles as collinear or usoft depending on whether p ^ 0. 
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and expanding about 5. 
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em 4 

where the IR divergence was regulated by dimensional regularization. To reproduce this 
result from the SCET computation requires adding contributions from regions A, B, and C. 
By using Eq. ([U|) we can add these contributions and still integrate over the full phase space 
in x 2 . For each region we can also expand the prefactor in Eq. (|73[) in the appropriate manner 
without effecting the LO results. Since x\ is fixed we need not worry about subtractions 
involving this variable. Overlap occurs from |Ab| 2 with region C if x 2 = and with region 
A if X2 = 1, so we find two zero-bin subtractions for these contributions. We multiply by 
(1 — Xi) 1+e to give the same normalization as Eq. (|76|). The regions with a soft particle are 
unsubtracted and give 
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1 - e 1 



K c = (l-a;i) 1+e / dx 2 - — -. — — = - + 0(5) . (77) 
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For region B with three collinear particles the naive contribution and its two subtractions 
(x 2 = 1, x 2 = 0) are 
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For the collinear integral we therefore find 

K B = K B - K B - K 2 B = ~ + 0(5) . (79) 
This result is IR finite as expected from the fact that the IR divergence comes from the soft 



region A in Fig. |TT| and not from the collinear region B. 
Adding the contributions from the three regions we find 

K A + K B + K C = -—-- (80) 

em 4 

in agreement with the result for K at this order. Thus with the zero-bin subtractions the 
sum of SCET diagrams reproduces the expected result for the 1/(1 — xi) bremsstrahlung 
term in the cross section da/dx\ as X\ — > 1. 
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FIG. 12: Contraction of fields in the time ordered product of LO SCET currents, T{J^\ J^}. 

Finally, we briefly remark as to whether we could have turned the phase space computa- 
tion into the imaginary part of a loop graph that we already know how to deal with from 
zero-bin examples in previous sections. One might think that the total cross section can be 
obtained by computing a forward scattering loop diagram in SCET and taking the imaginary 
part. However in some cases the optical theorem must be applied with care due to the mo- 
mentum scaling of different types of SCET fields. A simple example is the Born cross section 
for 7* — ► qq which we can consider computing from lm[i Jd 4 x exp(— iq-x)(0\ J(0) Jr(a;)|0)]. 
In the full theory there is a contribution with hard loop momentum and the imaginary part 
contributes to the total-u. In the SCET we are focusing on corners of phase space like 
back-to-back jets in the n and n directions. The product of LO currents allows for a 



loop with usoft momentum, shown in Fig. |l2|a, without violating momentum conservation. 



However this graph evaluates to zero due to the multipole expansion on collinear lines, 

d D k 1 

(2tt) d {n-{k + p)){n-{k + p) ( ' 

The product of currents does give a nonzero contribution, just not from the imaginary part of 
this usoft loop in SCET. Instead, the SCET fields give the imaginary part of two propagators 



as in Fig. 12b, yielding an integrand for the phase space integral that is accurate in the desired 



phase space region for each line. Because the matrix element factorizes into a product of two 



matrix elements, Fig. |12| b is not a disconnected contribution that can be discarded. This 
reduces the problem back to squaring the current, as depicted in Fig. [L?|c. It also gives a 
hint as to why the derivation of factorization formulas from SCET is more predictive than 
requiring a strict OPE in QCD, much as for diagrammatic factorization |38|, |72|, [73], [7JJ. 

VII. ZERO-BIN SUBTRACTIONS IN SCET n (EXCLUSIVE PROCESSES) 

In this section, we consider an SCET with degrees of freedom which are suitable for 
describing exclusive QCD processes with both energetic and soft hadrons. This theory is 
usually called SCETn, and contains fields that describe nonperturbative collinear and soft 
momenta as pictured in the p + -p~ plane shown in Fig. O. 16 It is also necessary to include 



16 



We do not need messenger or soft-collinear modes [[75f for the reasons discussed in section Vll A 
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FIG. 13: Degrees of freedom and momentum regions for SCETn which describe non-perturbative 
fluctuations with p 2 ~ Aq CD . The modes include n-collinear (c n ), soft (s), and when applicable 
also n-collinear (cn). In perturbation theory these modes extend all the way in to zero-momenta. 
For QCD the solid (red) curve represents the region where infrared divergences are rendered finite 
by Aqcd • Also shown (in pink) are three regions of perturbative momenta, two with hard-collinear 
momenta (hc n , hen) and one where the momenta are hard. 



these momentum regions when considering mixed inclusive and exclusive processes. We 
begin with a discussion of the ways in which SCET n differs from the SCETj and NRQCD 
examples discussed previously. 

In perturbation theory with massless particles, physical IR divergences occur as p 2 — > 
either with collinear scaling {p + ,p~) ~ Q(i] 2 ,l) or soft scaling {p + ,p~) ~ QiViV) f° r 
small dimensionless power counting parameters 77 and a large momentum scale Q. This 
is well known from the study of the Landau equations and use of the Coleman-Norton 
theorem [72, TfJ. In QCD, IR divergences either cancel between diagrams or are cutoff by 
the nonperturbative effects that generate confinement at a scale p 2 ~ Aq CD . In Fig. [13] we 
show the confinement scale by a red solid line. To formulate the power counting for these 
non-perturbative momenta, we take rj ~ Aqcd/Q- The collinear and soft fields represent 
distinct IR sectors as given by their momentum scaling, and together cover all approaches to 
the solid (red) curve in the p + -p~ plane. 17 The different sectors are separated by perturbative 
rapidity gaps. 



17 In this paper we do not consider processes which have important contributions from potential momenta 
for forward n-n scattering, k + k~ <C k\, which are sometimes referred to as Glauber modes. 
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The distinction between soft and collinear regions can be made with the variable 



C P = ^ , (82) 



which provides a measure along the solid red curve in Fig. |13|. For the different regions we 
have: 



n-collinear 
soft 

n-collinear 



C P ~ » 1 , (83) 



Thus to avoid double counting, we must make sure that in the variable ( p the n-collinear 
mode does not double count the soft mode and vice-versa, and also that the n-collinear 
mode does not double count the soft mode and vice-versa, etc. Note that the variable ( p 
provides a way of distinguishing the modes and at the same time allows us to maintain the 
boost-inversion symmetry |75|, 77|. The boost-inversion symmetry allows one to swap the 



soft and n-collinear fields, etc. when setting up the modes for the description of a physical 
process. 

As discussed in Ref. [|J, momentum conservation strongly constrains the form of soft- 
collinear interactions in SCETn- Adding a soft mode p s ~ Q{ViV>v) to a collinear mode 
Pc ~ Q{v 2 y 1) v)i which both have p 2 ~ Q 2 i] 2 , gives an offshell hard-collinear momentum 
Phc = Ps + Pc with p 2 hc ~ Q 2 i]. Thus all physical interaction Lagrangians and operators in 
SCETn will have > 2 soft fields and > 2 collinear fields. As long as double counting (and 
divergent convolutions) are avoided, we can group like fields together in gauge invariant 
products to obtain factorized amplitudes at any order in the power expansion in rj. In some 
cases, one can more directly prove that the convolution integrals converge J75| , and for these 
cases it is less important to be careful about the zero-bins. 

We will show that avoiding double counting in SCETn involves zero-bin subtractions 
similar to the previous sections, with the added complication associated with ensuring that 
regions in ( p are treated correctly. Due to UV divergences in rapidity, this requires a reg- 
ularization method. It also requires extra renormalization parameters for the insertion of 
any operator that connects soft and collinear fields, which we denote by « + and /i_. 18 In 
dimensional regularization, the parts of the SCETn action that are purely soft or purely 
collinear have the standard /i 2e multiplying couplings, so all factors of a s are a s («). The 
factors of /i + and /i_ only occur from mixed soft-collinear operators. We will show below 
how fi± appear in dimensional regularization, and also with a cutoff regulator. The fi + and 
/!_ parameters are tied together with the usual \x by the dynamics of factorization, which, 



Is 



For typographical convenience we use and /i+ although superscripts would be more appropriate. 
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independent of the UV and IR regulators, gives 



jjL = (JL + fi- . (84) 



Under an RPI-III transformation on the basis vectors, n — > e a n and n 
longitudinal boost on coordinates and fields), /i + behaves like a momentum p + = n-p and 
/i_ behaves like p~ = n ■ p. Furthermore ( p scales under a RPI-III transformation. These 
boosts correspond to a universal shift of all degrees of freedom along the solid red curve in 



Fig. [13], and thus do not change the fact that having distinguished between modes using ( p 
in one frame we also avoid double counting in any other frame. For a process with only soft 
and n-collinear modes, the boost-inversion symmetry allows us to interchange the role of 
these modes [[HJ. In Fig. [13], we boost to lower p~ and increase p + , so that the c n overlaps 
the s, and the s overlaps the c n . We then switch our definition of plus and minus, p + «-> p~~, 
with the overall outcome that c n <->■ s. Differentiating between modes using the variable ( p 
keeps them distinct throughout this process. 

The basic structure that we have in mind for a factorization formula in SCETn is 

J dk + dk~ dp~ dp + J(k ± ,p ± , /i±, /i' ± ) 0„(p", /i 2 ) s (/c + , k~ , //+, //_) (pn(p + , £t' + , /U 2 ), (85) 

where J contains perturbative contributions from both hard-collinear and hard momenta 
(as shown by the solid pink dots in Fig. [13]). For cases where only the n-collinear and soft 
modes are relevant, we have the slightly simpler form 

Jdk + dp~ J(/c + ,p~,/i+,/i_) n (p~,/i_,/i 2 ) (p s (k + ,n + ,fx 2 ) , (86) 

with J purely hard-collinear. If only the n-collinear and n-collinear modes are relevant, we 
have 

Jdk + dp~ J(p + ,p",/i_,/i' + ) n (p~,/i_,/i 2 ) 0n(p + ,// + ,^ 2 ) , (87) 

with J having hard momenta. The idea is that due to the separation of degrees of freedom 
in rapidity space the distributions can depend on /i ± . This dependence is similar to that 
for fragmentation functions in Ref. ]33|, |34j]. The meaning of the fi + and /z_ variables in 
the distribution functions is described further below in section |V11 (J| below Eq. (129). The 



presence of the fi± parameters allows us to formulate the non-perturbative matrix elements 



that give <j) n and S as boost invariant objects. This evades an argument made in Ref. [77 
that no IR regulator will allow a boost invariant factorization of soft and collinear modes in 
SCETn . Our proposed factorization formula differs from the conclusion of non-factorization 



in Refs. [[fvj |7l|, pOfl - The effects due to and fi + are actually not IR sensitive: they 
denote a choice we have to distinguish the IR regions. They behave like the dimensional 
regularization parameter /z in that we can compute the dependence on these parameters in 
perturbation theory because of the large rapidity gaps. 
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The soft and collinear modes in SCETn generate the physical hadron states in the effective 
theory, with each mode generating the physical states in its sector. The two sectors are 
separated by a perturbatively large rapidity gap, so we do not need to consider hadrons 
made of both modes. Since near the mass-shell p + ~ P 2 ±/p~, the rapidity scaling in Eq. fl85p 
gives gaps in p±/E of spacing Aq C d/Q- If we try to generate a hadron with an interpolating 
field built from soft and collinear fields, such as one collinear antiquark and one soft quark, 
then there are no physical non-perturbative poles by momentum conservation. Thus, the 
Hilbert space of states in the soft and collinear sectors are individually complete. 

In section |VII A| we begin by discussing a loop integral in SCETn taking into account the 
zero-bins. Our first one-loop example uses a hard cutoff regulator, in subsection VII A 1 



In section [VII B| we formulate the separation of soft and collinear modes using dimensional 
regularization, and repeat the one-loop example in subsection [VII B 1| . In section [VII C| we 



give a general discussion on how the zero-bin subtractions work on singular hard kernels to 
give what we call 0-distributions (the complete definition can be found in this section). In 
section [VII D| we apply this formalism to obtain a result for the 7*p — > it form factor at 
large Q 2 which is free from convolution endpoint singularities. In section V11E we discuss 
the so-called "soft" form factor for B — > n transitions, ( Bn (E), to argue that SCETn yields 
a result in terms of individual B and tt distribution functions. 



A. A SCETn Loop Integral with Subtractions 

As our first SCETn example, we consider a one-loop integral for the process U B — » ~f£i?" 
with Ery ^> Aqcd but using scalar quarks and gluons. The LO factorization formula for the 



full QCD process was considered in Refs. [^T], |82|, |84| using SCET. It involves n-hard- 



collinear fields and soft fields, but does not suffer from the subtleties in SCETn we wish to 
address. The toy example with scalars was considered in Ref . , where it was pointed out 
that this process with scalar quarks does not factor (naively) into a product of scalar and 
collinear terms, due to endpoint divergences which connect the soft and collinear matrix 
elements. This issue only shows up at subleading order for fermions. 

In Ref. J77], ^(J it was independently concluded that the convolution divergences en- 
countered in these situations spoil factorization. The analysis was based on a different IR 
regulator, implemented with an offshellness, and adding to SCETn a so-called messenger 
or soft-collinear IR regulator mode which has p 2 ~ Aq CD /Q. Any long distance colored 
interaction in SCETn would violate confinement in QCD and therefore be forbidden, but 
in perturbation theory one is free to introduce modes with p 2 below Aq CD if they facilitate 
the regulation of IR divergences. 19 The fact that messenger modes should be considered as 



19 In particular, any mode that would leave both soft and collinear modes onshell must have p 2 <C Aq CD by 
a parametric amount J75[ . In QCD a physical IR cutoff is provided by confinement which eliminates any 
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part of the IR regulator was discussed in Ref. |65], where it was shown that they are absent 



with an energy dependent gluon mass IR regulator, but that a common dependence on this 
regulator still appears in the soft and collinear matrix elements. 20 Messenger modes were 
also absent with the calculations using analytic regulators in Ref. |78[ |. 

With our definition of modes in SCETn in Fig. [13], messenger modes are not needed. If, 
on the other hand, one were to take the picture for SCETi in Fig. ^a and translate it into 
modes for SCETn in a one-to-one correspondence, then the usoft mode in SCETi becomes 
a messenger mode for SCETn f75|| . So the SCETi mode decomposition seems to want a 
messenger mode in SCETn. The fault here is with the translation, which required both a 
boost and a scale transformation. There is nothing wrong with the boost, but QCD is not 
scale invariant. Due to the presence of Aqcd ; translating the IR tiling of modes in SCETi 
into a tiling of IR modes for SCETn, gives modes which hide the physical situation, in 
particular the perturbative split of the hadronic physics in rapidity space. In our definition 
of SCETn the IR regions that show up in perturbation theory and were described by the 



messenger mode in Ref. |75[ are absorbed into the soft and collinear fields. 

In considering scalar ll B —* / ~f£u 11 we are really treating a subleading contribution to the 
physical process with fermions where the photon comes from fragmentation of qq — * 7 or 
from a subleading contribution to the direct 7-production. In the fragmentation case we 
will have a non-perturbative soft distribution associated with the initial state B, <p s (k + ), 
and a non-perturbative collinear distribution for the 7 associated with the fragmentation, 
^n(p )- To factorize this physical process we must consider the imaginary part of the forward 
scattering, since due to the probabilistic interpretation for the fragmentation function, there 
will be a factorization formula for the decay rate but not for the amplitude. Examples of 



one loop diagrams in the full theory and SCETn are shown in Fig. |14[ The SCETn graphs 



in Fig. |14| b correspond to direct production, while those in Fig. |14| c are fragmentation. The 
SCETn factorization allows these two effects to be distinguished. 

In the full scalar theory we have a charged scalar "b-quark" with field (p b , charged light 
scalar "u-quark", 0, and neutral "gluon" and "photon" fields <fi 9 and (p respectively. The 
"weak" current and interaction terms are 

JZt = , 4uli = + , (88) 



such modes |8f| as they would physically correspond to colored degrees of freedom propagating between 
color singlet bound states that have already hadronized. They have been termed hyper-confining modes 
by Rothstein |8(J. The soft-collinear messenger mode considered in Refs. j7j| [77J [8(J are in the hyper- 
confining category. 

20 Note that in general matching computations between QCD and SCETn can be a bit tricky because one 
must be sure that ones choice of IR regulator is treating the Hilbert space of full QCD in exactly the 
same way as the Hilbert space of SCETn • An example of this type that we encountered in SCETi was 
our discussion of the field redefined current J'( ' in b — > sj in sections VIA and VI B. For SCETn further 
discussion of this point is given in Appendix |b]. 
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a) full theory graphs 





b) SCET n soft graphs 
P 





c) SCETjj collinear graphs 



P 



P 



-k -k 

FIG. 14: One loop cut-graphs for the 7-fragmentation contribution to B — > 7^2/ in the full and 
effective theories. The ®'s denote an insertion of the weak current and the leptons are not shown. 
The graphs in a) are in full QCD, while the graphs in b) and c) are in SCETn • For b) we have 
soft propagators, while for c) the propagators are collinear (dashed). 



with standard kinetic terms for ^ with mass m&, and for the massless charged and neutral 
scalars. Here the coupling G tracks the current, and g and e are coupling constants of mass 
dimension one. In SCETn the leading order Lagrangian is split into soft and collinear fields, 



C 



(0) 



2m b <j% iv-dcjil + cfil {id) 2 S + - 0f {id) 2 0f + g [0* <f>\ t <f> a>k ] 



(89) 



+ 0t {td) 2 <Pn+\<PZ (^) 2 <f>l+±<% {td) 2 <% + g cf>l p <j> n , q ] + e [<t>l P _ q 4 )P K, q \ ■ 



where cj) h v is a scalar HQET field [p8 |, the <p s and n are soft and collinear massless "quarks", 
and 4>9, 4>ni an d 4>Z are scalars for the "gauge" fields. The collinear fields have label mo- 
menta p" ~ rf and residual momenta pj- ~ 77 and p+ ~ r/ 2 , and the soft fields have label 
momenta k + ~ r] and residual momenta k^r ~ 77 and ~ 77. Writing out the label and 
residual terms explicitly in the kinetic terms we would have {id fl )<j) s — > {n^V /2 + id^)4> s ^, 
{id^cpn — ► {n^V '/2 + id^)(pn tP etc., with the standard treatment of leading and subleading 
terms. Demanding ~ 77° the power counting in 77 is 



3/2 



77 . 



V: 



<t> 9 s ~ 4i ~ V • 



(90) 
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The leading order currents we will need are 

j(0o) 

[n-pn-l— iO 



t(0o) 



n-pn-£—iO + )(n-pn-k—iO + ] 



op>= ^ 



j(0c) 

ff C) = J] 7^ 5 - n+ w- „ - n+ N [<rf,-^ 0« 0n,p-g 0n -J , 

where np ~ 77 while r*£ ~ r}k ~ 77. Note that the operators and include restrictions 
on the label sums with n-i 7^ and n-p 7^ respectively. In these bins the l/(n-pn-£) factor 
would be a collinear or soft propagator in SCETn, and these bins are taken into account by 
time-ordered products with subleading SCETn Lagrangians as discussed in Ref. p7| . 

Note that in scalar SCETn, we can construct leading order operators with additional 
scalar fields since the extra powers of 77 are compensated by 1 jin ■ p n ■ t) factors and the 
mass dimension is compensated by the dimension of the couplings. In our perturbative ex- 
ample only the currents shown are needed (plus counterterm operators). A non-perturbative 
treatment would require additional terms. In the scalar theory, matching tree level graphs 
gives = eG, = egG, and J ( - 0c - ) = g 2 G. In dimensional regularization, the currents 
are modified in a manner described in section IV11B. 



In the gauge theory with fermions, the full set of operators for B — > yev would be 
determined by gauge invariance and leading order matching. The current analogous to 
is LO, while the currents analogous to ojj 6,0 ^ are suppressed by one power of 77. 
The non-perturbative treatment is simpler in gauge theories since it is constrained by more 
symmetries. See for example, the construction of soft-collinear SCETn operators in Refs. 173, 



A goal of this section is to demonstrate that although there is not a simple factorization 
for the scalar B — > ■yev process of the form J(/i) ®<j) n (fi) <g> <j> s (p)i there appears to be a more 
involved factorization which contains terms of the form 

Jdk + dp~ J(k + ,p~,n + ,n_) <fr n (p~,n_,fi 2 ) s (fc + ,/i+,/i 2 ) . (92) 

Here J is a perturbative jet function, and the 0's are non-perturbative, with (j) n given by a 
matrix element of collinear fields and <ft s given by a matrix element of soft fields. In Eq. (|9"2"D 
we have two additional scale parameters /z + and /i_ in SCETn (or two factorization scales 
in a more traditional language in full QCD). We will demonstrate that these scales are 
connected in a specific way to the standard renormalization scale 

H- H+ = ^ ■ ( 93 ) 
Although our conjecture about the existence of a factorization formula valid to all orders in 



a s differs from Refs. [77, 78, 80], the structure of the result is also different from standard 
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factorization formulas in the literature. Because we do not have the identity /i + = /i_ = //, 
there is not a simple factorization for the amplitude for this process, in agreement with 



the conclusions in Refs. J77, |7J| 5D| . The perturbative formula that we find for the observ- 
able process, Eq. (|85|), has all the desired properties of a factorization formula, including 
correctly reproducing IR divergences (from the zero-bin subtractions), finite convolution in- 
tegrals (from the zero-bin and renormalization), and distinct matrix elements for the soft 
and collinear objects. 

In the context of the analytic IR regulator used in [78], it was pointed out that there 
was an interesting cancellation between IR divergences in soft and collinear diagrams. Our 
observation is that this cancellation has to do with avoiding double counting just like the 
zero-bin subtractions, rather than having to do with reproducing IR divergences in QCD. 
For physical observables, like the forward scattering graphs, the method for avoiding the 
double counting is computable, and can be handled in perturbation theory. It results in 
regularization parameters \x+ and which encode the coupling between soft and collinear 
modes with a simple correlation. Our results turn situations which were previously plagued 
by the unphysical convolution endpoint singularities, into manageable finite amplitudes, 
which one can then try to arrange into a predictive factorization formula. Since the divergent 
effects are computable they do not spoil many of the nice features obtained in simpler QCD 
factorization formulas. 



I. Soft- Collinear Division with a Hard Cutoff Regulator 

We begin by considering a hard cutoff, a, between the soft and collinear modes, as 
indicated in Fig. [U| For a loop momentum fc M we define ( k = n ■ k/n ■ k = k~/k + as 
discussed near Eq. (|T5|). Only the magnitude of the rapidity variable, |Cfe|, is relevant for 
distinguishing the soft and collinear modes. Switching variables from {k + ,k~} to {k + ,(k} 
gives dk~ = \k + \d( k when we integrate over — oo < k + < oo and — oo < ( k < oo, so the loop 
integral is no longer analytic in k + but remains analytic in and likewise if we switch to 
{k~Xk}- When imposing hard cutoffs we need to avoid the physical poles, which can be 
accomplished using cutoffs in Euclidean space after Wick rotation. For variables {k + ,(k} 
the Wick rotation k~ — > ik~ is equivalent to ( k = K'k-> while for {k~,(k} the Wick rotation 
k + — ► ik + gives ( k = — i(' k . In our examples Wick rotation about the origin suffices, and 
the poles in complex (' k occur along the imaginary axis in the first and third quadrants. We 
take cutoff's 

soft: - a 2 < C' k < a 2 , (94) 
collinear: — a 2 > ( k or (' k > a 2 . 

As mentioned above we only require n-collinear and soft fields in SCETn for the example in 
this section and so are free to include the entire Cfc ~ 1 region in the soft modes. (For more 
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FIG. 15: Graphs with scalar propagators as a toy model for the 7-fragmentation contribution to 
B — > with fermions. The (g> denotes the weak current and the leptons are not shown. Graph 
a) is in full QCD, graph b) has an insertion of Oj? 6 ^ and a loop with soft fields in SCETn , and 
graph c) has an insertion of Oyf 1 and a collinear loop (with dashed propagators). 



complicated problems the region (' k ~ rf would need to be disentangled for the n-collinear 
modes.) We can take a 2 ~ rj, and note that under an RPI-III transformation on n and n |79[ 
(a longitudinal boost) that a 2 behaves like a {p~) 2 momentum just like (k does. Also note 
that the power counting scaling which fixes the soft and collinear components only depends 
on \( k \ or \(' k \ and so does not care about the Wick rotation. 

For simplicity we consider the same diagram as discussed in Ref. [78| which is shown in 
Fig. |T5|a. Unlike Ref. |[78| , we do not analyze this graph with the method of regions. Instead 
we consider the diagram in full scalar field theory (Fig. |l5|a) and the corresponding diagrams 
in scalar SCETn (Fig. |T5|b,c). The difference of the two results gives a matching contribution, 
and allows us to check that the full theory IR divergences are correctly reproduced with 
zero-bin subtractions implemented in SCETn . It also allows us to discuss the factorization 
formula in Eq. (]HB|). 

For simplicity we will leave off the prefactor ieg 2 G/(p~£ + ) in quoting results for graphs 
in this subsection. For the full theory diagram, we have the integral 

d D k 1 

(95) 



rscalar 
'full 



(2tt) d [(k - t) 2 + z0+] [k 2 + i0+] [(k -p) 2 + z0+ 
Evaluating this onshell with p M = p~n' M /2, £^ = £ + n' M /2, so that p 2 = £ 2 - 



0, we have 



rscalar 
-"full 



16vr 2 (p-^ 



1 



-IR 



em 



P' 



+ iln 2 

2 



p~£ + 



7T 
12 



(96) 



Here the IR divergences are regulated by dimensional regularization. For the soft and 
collinear graphs in Fig. [I5|b,c we find 

d k r - 



rscalar 
'soft 



rscalar 



E 

fc+^O 

E 

fe-^O 



1 



(27r)^ [k 2 - 
d D ki 



-n-£n-k + i0+][k 2 +i0 + ][-n-pn-k + i0+] 
1 



(97) 



(2tt) d [-n-£n-k + iO+}[k 2 +iO+][k 2 -n-pn-k + i0+] 



where in both cases the first two terms are the displayed propagators, and the last factor 
comes from the non-local vertex which emits the scalar soft or collinear fields in SCETn. 
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To compute the EFT graphs we implement the hard cutoff in Eq. ( |94]) to regulate UV 
effects in the effective theory diagrams. With this regulator the zero-bin subtractions are 
automatically zero since they are outside the region of integration. The hard cutoffs are 
theta functions in the integrand so they give identically zero for the integrand evaluated 
in the subtraction regions. Therefore with this regulator the full integrals are given by the 
naive replacement in Eq. fliCf). We discuss in detail the calculation of the SCETn diagrams 
in Appendix |A]. For the soft graph the result is 



rscalar 
'soft 
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lQn 2 (p~£+) 
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em 
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(98) 



where we defined fi + = fi/a. Note that since a boosts like a minus-momentum, /i + behaves 
like a plus-momentum, and the result in Eq. (p8|) is RPI-III invariant. For the collinear 
graph the result is 

-i r i 

— — In ( — )-!- Iir | 
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" 16vr 2 (p-^ 
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16 
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(99) 



where we defined /i_ = a/i. Here /i_ behaves like a minus-momentum and the result in 
Eq. (^) is also RPI-III invariant. The soft and collinear regularization parameters fi± 
defined in the computation of Eqs. (|98|) and fl99|) obey the anticipated relation, 

■ 2 (100) 



where the a dependence cancels out in this product. Moreover, with a 2 ~ rj we find that 
with /i 2 at the matching scale, /i 2 ~ QAqcd, one can still take fi + ^ Aqcd and /i_ ~ Q. 
Thus we can simultaneously minimize the logarithms in the SCETn matrix elements, which 
in our perturbative computation are represented by I^ 3,1 and J^ alar . 

rscalar + L s ^ alar , We find 



Adding the two SCET n graphs, I ; 
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16n 2 (p-£+) 

We see that with the relation ii 2 = yU+/i_, the 1/eiR poles agree exactly with the full theory 
expression in Eq. (|96|) as required. To match the full and effective calculations we set 
" " and subtract to find 
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(102) 
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Here the ln 2 (p - /! + /£ + /i_) contributes p~ and i + dependence to the jet function 
J(£ + ,p~ , //+, at one loop. To minimize the large logarithms in the matching calcu- 
lation in Eq. (|102|) we take /i 2 = to be of order the hard-collinear scale, and take 
~ rj. Since the matching result in Eq. (|102| ) depends on yU + //i- it depends on a, 
which is not surprising. Here a ensures there is no double counting between the soft and 
collinear modes in the IR, but a also changes the behavior of the collinear and soft modes 
in the ultraviolet. This change is compensated by the perturbative Wilson coefficient, and 
in perturbation theory the sum of these contributions reproduce the full theory result. 

The result in Eq. ( pL02| ) is shown for illustration only, since a complete matching calcula- 
tion for scalar B — > 'yiv requires a computation of all diagrams, not just the one diagram 
that we considered. For example, one should also compute graphs with the scalar gluon 
attached to the 6-quark line, and wavefunction renormalization type diagrams in both the 
full and effective theories. 



B. Dimensional Regularization Division in SCETn: General Discussion 



In this section, we discuss the use of dimensional regularization for the UV divergences and 
the separation of soft and collinear modes. This regulator makes higher order computations 
more feasible and preserves gauge symmetry. We also expect that it will make it easier to 
compute anomalous dimensions and sum logarithms using renormalization group techniques, 
although we do not address these features here. Finally, it is useful to consider dimensional 
regularization in order to compare how the separation of modes in rapidity space appears 
with a different regulator. 

Since the standard application of dimensional regularization is boost invariant, it does not 
provide the ability to distinguish modes in rapidity space. This also means that in general, 
divergences in the rapidity will not be regulated by standard dimensional regularization. 
For an insertion of a mixed soft-collinear operator, we can regulate the rapidity space in 
dimensional regularization by scaling out factors of the label operators from the Wilson 
coefficients. To implement Fig. [13| in dimensional regularization, the correct form of the 
operators are 



J(pJ, k+) [(q s S) kt T s (Sks) kt ] [{inW) p -T n {W^ n ) p -] 



(103) 
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Here the label operator V gives the plus momentum from soft fields, and the label operator V 
gives the minus momentum from collinear fields. The momenta subscripts occur for products 
of quark fields and Wilson lines, (S^q s ) k+ = 5(k + -V)(&q a ), (W^ n ) p - = 6{p~ -P)(W%), 
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which ensures that the momenta are gauge invariant and that the gauge symmetry is not 
spoiled by the factors of |P| e , etc. The absolute values ensure that we raise a positive 
physical momentum to the e power, and thus do not modify the cut structure of matrix 
elements. 21 J is the Wilson coefficient jet function. 

This rescaling will allow us to properly distinguish the soft and collinear modes in dimen- 
sional regularization without imposing a hard cutoff to implement the division in Fig. [13. 
This modification of the current is not done to solve a problem in the IR - it is the zero-bin 
subtractions for the soft and collinear fields which will ensure that there is no IR double 
counting. The zero-bin subtraction terms are integrated over all space, which introduces 
new UV divergences in rapidity space, and in Eq. (|103|) the factors of \V\ e etc. are necessary 
to regulate these UV divergences. If one thinks of splitting the full loop integral / into a 
naive part / and a subtraction part Iq, then I has an IR rapidity divergence, while Iq has 
both UV and IR divergences. The IR divergences cancel in I — I — Iq, so the rapidity di- 
vergence in / is pure UV. We will see that these remaining UV divergences can be removed 
by counterterms. 

Before giving the rules for constructing Eq. ( |103| ), let us consider how it should be used. 
When we do a collinear loop involving an insertion of this operator we expand in 
but DO NOT expand the {k + / factors, and we do the opposite for a soft loop. This 
dimensional regularization rule is forced on us in any field theory with a multipole expan- 
sion, and SCETn has a multipole expansion between components of the soft and collinear 



momenta. The rule was discussed in Ref. p3] for NRQCD in examples involving mixed 



usoft-soft loops. In general one does not expand matrix elements of the soft fields in D- 
dimensions when doing the collinear loops and one does not expand matrix elements of 
the collinear fields when doing soft loops. The factors of (p~//i_) e and {k + / fi+) e should 
be thought of as being associated with the renormalized coupling function J, just like a 
factor of /i e is associated to the strong coupling #(//). For purely collinear or purely soft 
operators we apply dimensional regularization in the usual manner. The only place that n± 
appear is in the insertion of a mixed soft-collinear operator. All purely soft operators and 
purely collinear operators only have fi e factors, and so all couplings are a s (fi). The couplings 
do not dependent explicitly on /i ± . In multiloop diagrams one can carry out the standard 
renormalization procedure first, and leave to the end the rapidity renormalization for the 
final loop involving the soft-collinear vertex. 



Recall that the labels are positive for particles, and negative for antiparticles fcs]. Combining both particles 
and antiparticles into a single field distinguished by the sign of the label simplifies the formulation of the 
effective theory. One could instead have used separate fields for the particles and antiparticles, in which 
case the antiparticle field could also be chosen to have a positive label. The absolute values in \P\ e mean 
that we are using the momentum of the particle, which is unambiguous, rather than the label on the field, 
which is convention dependent. Due to the zero-bin conditions 7^ and 7^ 0, there is no problem 
at the origin. 



54 



Lets consider how we determined the powers of //+, //_ and \x for Eq. (|103|) . In dimensional 
regularization, factors of fi e appear from ensuring that coupling constants in the action are 
dimensionless. Demanding that this is the case for J gives the /z 2e factor. To determine 
the factors of [i± we must examine the scaling of fields in the operator along the solid red 
curve in Fig. |T3|. An RPI-III transformation scales all modes by a common amount and 
separate invariance under this transformation demands we introduce a to compensate 
each V, and a for each V. Charge conjugation requires the same power for quarks 

and antiquarks, and the boost-inversion symmetry [75, |77J requires the same power be 
used for the soft and collinear fields. We can also demand rapidity invariance under small 
individual scalings of the soft and collinear sectors. This will determine the power of the 
label parameters, \V\ e and \V\ e . The power is related to the space-time dimension because 
in doing this rapidity scaling we demand that the invariant mass p 2 remains homogeneous. 
As an example consider a scaling by (3 > 0: Xn, P -(0) — > Xn,Pp-{fy = P~ t Xn, P -(fy, where 
here x = as in the soft-collinear operator. To derive the (3~ e factor write the quark field 
Xn , Pp -(0) = 5(j3p~ - P)Xn(0) = J d d k 5(f3p--k-)5(k 2 )6(k~)a k , and then shift k~ -> (3k~ 
and k\ — > (3k"j_ (to keep k 2 homogeneous). This f3~ e factor from the transformation of Xn,p- 
is exactly canceled by the \V\ e — > \V\ e /3 e factor acting on this field. The symmetry of the 
problem dictates that we need one such factor for each collinear field. Repeating these 
arguments with a scaling parameter (3' in the soft-sector determines the \V\ e terms. 

These arguments determine the proper operator for the dimensional regularization com- 
putations in SCETn, with an example shown in Eq. ( |103[ ). The Wilson coefficient J has 
non-trivial //_ and /x + dependence which cancels the dependence on these parameters in the 
matrix element of the effective theory operator order by order in ot s {p). Thus we see that 
factorizing the soft and collinear modes in SCETn also requires introducing fi + and just 
like with our cutoff regulator. 

It is interesting to compare the regulator introduced in Eq. ( p.03p with the use of analytic 



regulators used for collinear computations in Refs. (7^, [89], 90]. Much like an analytic reg- 
ulator, the result in Eq. ( |103|) modifies the power of a momentum dependent factor in the 
integrand. It is used to regulate divergences that are not handled by dimensional regular- 
ization, which is also the motivation for introducing an analytic regulator. However, unlike 
the use of analytic regulators, Eq. ( |103|) is gauge invariant, defines the modification at the 
operator level, and does not modify the power of the propagators in the EFT. Furthermore, 
as already emphasized, with our zero-bin setup this regulator is needed for divergences in the 
UV, rather than the IR. Since these divergences arise due to the separation of momentum 
fractions in hard scattering kernels and collinear operators, we anticipate that the addi- 
tion of a power of the momentum fraction for each labeled field will regulate UV rapidity 
divergences in a general situation. 
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1. Dimensional Regularization for the one-loop example 



We now repeat the computation in the last section with dimensional regularization for 
the UV. The loop integral in Eq. (|95|) has IR divergences for k — > 0, k — > £ + , and k — > p~ . 
We take one of the propagator lines to have an infinitesimal mass m 2 to regulate these IR 
divergences. For simplicity we will leave off the prefactor ieg 2 G/(p~£ + ) in quoting results. 
For the full theory diagram we have 
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(104) 



Here m 2 regulates the IR divergences in a manner similar to the solid red curve in Fig. [13|. 
Other choices of IR regulator can be made, and in Appendix || we repeat the computations 
done in this section with i) factors m|, m 2 , and m\ in the three propagators in Eq. (|104|) , 
and ii) taking p 2 ^ and £ 2 ^ in Eq. (|104j) . The choice mi = is also discussed in 
Appendix [B], but makes the matching more complicated. 

The LO mixed soft-collinear SCETn currents in dimensional regularization include the 
UV rapidity regulation factors, and are 
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where we have suppressed the sums over label momenta shown in Eq. fl9lD , and in general 



the J® are functions of the label momenta (p , £ + , k + , . . .). Using the currents Ofj and 
0{j°' ) for the soft and collinear graphs in Figs. [i~5|b,c respectively we have 
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Here the k + ^ and k ^ conditions denote the overlap regions where the soft integration 



variable becomes collinear and the collinear integration variable becomes soft, as in Fig. [13. 
The sums over k + ^ and k~ ^ ensure that the [— p~k + ] and [— £ + k~] propagators never 
get small. By examining the scaling, we find that no subtraction is necessary for k + ^ £ + 
and k~ ^ p~ here, so though present, these restrictions were not shown. Eq. (|TTD tells 
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us that unlike the SCETi computations and the SCETn cutoff computation, here we have 
zero-bin subtractions for both the soft and collinear diagrams. These will ensure that we 
do not get spurious singularities from the [— p~k + ] and [— £ + k~] propagators. The naive 
integrals and subtraction integrals are 
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Note that we must keep the m 2 dependence in the subtraction integrals to properly avoid 
double counting the zero-bin regions in the differences — ^osoft r anc ^ IIT X&V ~ ^ocn la ' " i 

which from Eq. flTTD give the result for I^ft^ anc ^ IqT 1 ™ respectively. 

For the soft graph we do the k~ integral by contours. Due to the pole structure this 
restricts the fc + -integration to the region < k + < £ + . The k± integral is then done. 
For the soft subtraction integral we follow the same procedure which this time leaves the 
integration region < k + < oo. We find 
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In the last line the first {■ • • } factor comes from the integral over < k + < £ + , and the second 
from £ + < k + < oo. Computing the full soft integral in Eq. ( |106| ), /foft ar 
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Much like the examples in SCETi the zero-bin subtraction integral /osoft 1 cancels the IR 
singularity in the k + integration in I^f^ and replaces it by a UV divergence. 
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For the collinear integrals, we do the contour integration in k + which restricts the re- 
maining integration region in k~ . For the naive and subtraction integrals we find 
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which for our example, are the same integrals as for the soft loops but with £ + 



find 



Thus for the complete collinear result in Eq. ( |106[ ), I\ 
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The results in Eqs. ( |109|) and (|111|) have l/euvhi(m 2 ) divergences, terms that did not 
appear in our example with a rapidity cutoff, and are simply artifacts of the dimensional 
regularization setup. These divergences arise from the fact that the UV collinear divergences 
induced by the zero-bin subtraction are multiplicative over all loops and propagators. They 
are canceled by a special type of counterterm that is proportional to the renormalized dis- 
tribution function at the origin, 0(0, /z)/euv- The presence of these operators is discussed 
further in section |V11 G below. Here we have one such counterterm current for the soft loop 
and one for the collinear loop 
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These correspond to counterterms for operators which give a 0(0, /i), corresponding to the 
initial B meson, and a 0(0, /a) for the quark part of the final photon wavefunction. Note 
that the limit k + —>■ and q~ —>■ is done at the end. 

The necessary counterterm coefficients for the results in Eqs. ( |109|) and fllllj) are 
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At one loop these operator generate graphs similar to the diagrams in Fig. |T5]p,c. Using 
the same IR mass regulator, pulling out the same prefactor as the other diagrams, and 
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performing the standard UV renormalization of the operators in MS prior to multiplying by 
the rapidity counterterm gives 
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Due to our choice of 5j(° d ' 0e ) these counterterm diagrams exactly cancel the 1 / 'euv hi(m 2 ) 
terms in the collinear and soft loops. In dimensional regularization, this type of countert- 
erm operator is quite important to the rapidity renormalization, as we discuss in the next 
section. In particular in going beyond perturbation theory, these same type of counterterms 
are required to cancel UV divergences in the convolution over the non-perturbative matrix 
element. 

Adding the soft, collinear, and I^od*' -^ctocT counterterms we find the SCET n result 
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The effective theory still has UV divergences shown on the second line of Eq. Q115 ). These 
divergences occur because of the separation of ( p momenta. The remaining UV divergences 
are canceled by a counterterm for the jet function coefficient j( 0a \ Putting back the pref- 
actor, we find the counterterm 
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As is familiar from SCETi, the counterterm depends on logs involving /x's. Thus, finally, the 
renormalized EFT result is 
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Comparing the first two terms in Eq. ( p.!7|) with the full theory result in Eq. (|104 ) we 
see that the IR divergences are exactly reproduced if and only if 



(118) 
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Thus again this condition follows from the dynamics. It is interesting to note that the 
\n(p~£ + ) ln(m 2 ) divergence is reproduced independent of the power of \V\ and \V\ used in 
Eq. ( |103| ), but that the ln 2 (m 2 ) term is only reproduced for the power e, which was derived 
in section |V 11 B[ Using /i 2 = the difference of the remaining finite terms gives a 

contribution to the one-loop matching 
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From Eq. ( |119[ ) we see that the jet function will be a non-trivial function of fi + and 
(and thus /i 2 = //+//-) whose //-dependences will cancel against dependence on these vari- 
ables in the SCETn matrix elements. The matching coefficient arises from integrating out 
perturbative effects associated with p~l + as well as perturbative effects responsible for the 
rapidity gap between the soft and collinear modes, and therefore has a different structure 
than what would be obtained if only the former were integrated out. In Appendix || we 
verify that the same result for /match ^ s obtained when we regulate the IR divergences in the 
full and effective theories with three non-equal masses, 7711,2,3, or when we keep the I 2 and 
p 2 offshell. The result should be the same because the matching only depends on the UV 
regulator which we keep the same in these computations. Note that in Eq. ( |102| ) we used a 
different UV regulator than in Eq. ( |119|) , which explains why the tt 2 terms differ. 

The results in Eq. ( |116| ) and (|119| ) are shown for illustration only, since of course, the com- 
plete anomalous dimension and matching calculations require a computation of all diagrams, 
not just the one diagram that we considered here for illustration. 



C. Zero-Bin Subtractions in Convolutions: General Discussion 



The remaining application of subtractions in SCETn will be for factorization formulas 
which appear to suffer from singular convolutions at the level of tree level matching. Much 
like the example discussed in SCETj, in SCETn we must avoid the zero-bin in hard scattering 
kernels defined by tree level matching. Doing so removes the double counting problem and 
renders singular convolutions finite. Here we only deal with the rapidity renormalization, so 
we make the simplifying assumption that the standard UV divergences have already been 
taken care of, and do not interfere with the steps carried out here. 

We will use dimensional regularization to separate the modes in rapidity space as in 
Eq. ( |103| ). To see why the convolution integrals are always finite, let us consider the vacuum 
to pion matrix element of a hard scattering kernel J(p~, fi~, fi 2 ) and a twist-2 collinear 
operator in SCETn that gives the light-cone distribution function (f) n (x,fi). For simplicity 
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we will not write the /z-dependences for J and (j) n below. 22 This leaves the matrix element 
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where we switched to dimensionless variables £i ;2 via p 1 = x\ n-p n and p 2 = x 2 n-p n , and in 
the second line we inserted the standard definition of the twist-2 distribution function 

<^(^)|w n , p -^75<_p-|0) = -if 7T 5(n-p 7T -p7-p 2 )<j) 7T (x 1 ,x 2 ,fi) , (121) 

where the 5 function gives conservation of momentum. Now suppose that we computed J at 
tree level (by a matching computation) and found that J(xi,x 2 ) = l/(p7) 2 = l/[(n-p n ) 2 x 2 ]. 
If we were not careful about the X\ 7^ condition, this would lead to a singular convolution 
integral as in Eq. (|2|). The zero-bin subtraction formula in Eq. (|17D tells us to impose the 
momentum conserving (^-functions carrying through all zero-bin constraints. Since the x 2 - 
integration is not singular, there are no zero-bin subtractions for x 2 7^ and we can combine 
the sum over label x 2 momenta and integral over residual x 2r momenta back into a integral 
over all x 2 momenta using Eq. ([IB]) : 
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where X\ — 1 — x%. If there had been zero-bin subtractions for x 2 they would carry through 
as additional zero-bin subtractions at x\ 7^ 1 after removing the ^-function. In the last 
line we set <fi n (xi,Xi) = 9(l — xi)9(xi) ^(xi) to make the support of the non-perturbative 
distribution function explicit. To turn the final sum over labels and integral over residual 
momenta into an integral over xi, there will be zero-bin subtractions from Eq. fllTf) . The 
subtraction acts on the integrand including the ^-functions, but just as in our perturbative 
analysis, it does not act on the |xi(l — xi)\ e factor. The expansion for x\ 7^ is from the 
right, about x\ = + , since this is how the variable scales towards the zero-bin region: 

„2 



9(l-x 1 )9(x x )^(x x ) = 9(x x ) [4(0) + Xl 4(0) + ^ 4(0) 

+ 9(x 1 )fa(0) [6(1- Xl ) +...]. 
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22 Note that integer powers of the p i can be moved from J to the collinear operator by inserting powers of 
V, but that our analysis is independent of this freedom. 
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In the set of terms obtained on the first line, the 9(1 — x\) disappears in the series so the 
support of the x% integration for the subtraction terms differs from that for the naive integral. 
This is the same as what we saw in our perturbation theory example in Eq. ( |110| ), where 
the naive integral was integrated over k~ G [0,p~], ie. X\ G [0,1], but the the subtraction 
integral was integrated over k~ G [0,oo], ie. x\ G [0,oo]. In the last line in Eq. (|123|), the 



terms are all zero (or finite subtractions) for the cases considered here, and therefore these 
terms do not contribute for our choice of zero-bin scheme as discussed in section [TV| 

Lets make the standard assumption for the twist-2 distribution that (j) n (0) = 0. Then 
using Eq. ([17]), the result for A w is 
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where as usual only the subtraction needed to remove the singular term was kept. Next we 
split the integration into a finite integral x\ G [0, 1] where the factor of |a:i(l — | e can be 
set to 1, and the integral of the subtraction term over x\ G [l,oo] where the e dependent 
term is needed 
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Here terms of 0(e) have been dropped. 

Eq. ( |125| ) is UV divergent, but we must still add to it the pion matrix element of the 
counterterm operator. This operator is determined by the UV counterterms that are neces- 
sary to renormalize our original operator, and can be derived in perturbation theory with 
any desired external states. Carrying out a one-loop computation with external quark states 
and using our perturbative kernel J = l/(p^) 2 we find the counterterm operator 
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with a counterterm coefficient 5C, 
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-l/(2euv)- The derivative with respect to (p^ +p 2 ) 
removes surface terms. In the vacuum to pion matrix element, they would result from a 
d/dpi of the 5-function in Eq. (|121|) if we had left out the d/d(pi +p%)- At tree level with 
quarks the matrix element of this operator vanishes — one obtains S'(p~) factors and the 
quark states have non-zero p~ momenta. The vacuum to pion matrix element of O^} gives 
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where the superscript (1,0) indicates a derivative with respect to the first argument. We 
should also include the matrix element of the finite part of the counterterm operator in 
Eq. (pq ) 0, which gives 

Af = -iA_c2V)&(0). (128) 



Adding the A c * 1 term to Eq. (|125 ), the UV divergence cancels, and sending e — > we obtain 
the finite result 

A n + A^ 1 — f^ M^^xMO,,) + 

n-p w { J [xiY \ /i_ 

= Cdx^t^- (129) 

n-Pn Jo (^1)0 



As indicated, performing the steps outlined from Eq. (|120|) to (|129|) defines the 0-distribution 



in dimensional regularization with our renormalization scheme. The //_ in the distribution, 
4>(xi, /i, /i_) is a short hand for the dependence on the ln(/i_) in the first line of Eq. (|129| ). 
Once again, in Eq. ( |129| ) the zero-bin subtraction has converted an IR divergence into a 
UV divergence - the naive IR divergence in the convolution has been converted into a UV 
divergence for the operator in Eq. ( |121| ), which is canceled by the operator renormalization 
counterterm in Eq. ( |126| ). Essentially the 0-distribution notation on a variable, (x)0 indicates 
that we have a sum over labels x 7^ 0, and do an integral over residuals dx r , together with 
applying the rapidity renormalization procedure outlined above for the UV divergences. The 
/i_ dependence in Eq. ( p9|) is canceled by Cj} ] (//_) ^ in Eq. (p8|). 

For other cases, the steps in determining the result for the 0-distribution are the same as 
in our example; however it should be clear that the final result will depend on how singular 
the perturbative kernel is, as well as the endpoint properties of the non-perturbative function 
that the 0-distribution is acting on. In particular, if the starting kernel was not singular there 
would be no zero-bin subtractions and we would obtain the naive result for the convolution 
that one finds without the 0-distribution. Note that if we had implemented a hard cutoff 



as in section |VII A 1| rather than dimensional regularization, then lower limits, like x > 5, 
would be induced on the convolution integrals, together with compensating 5 dependence in 
the jet functions. 

For illustration, we consider a few other cases in dimensional regularization that are 
quite common and which appear in the examples in the next section. First consider a 
distribution 0£(:ei, x 2 , that does not vanish at its endpoints, integrated against a kernel 
J = l/(n-p n p^). For the analog of Eq. ( |121[ ) we take the matrix element of the operator 
to give -if n fa8(n-p 7r - p\ - p 2 ) (f^(x 1 ,x 2 , fi) where fa = ml/ (m u + m d ), using Eq. (|142| 
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below with n ^ n. The steps leading up to Eq. ( |122| ) are very similar, giving 
.f n fi n (n-p^ 



B, 



—i- 



^2 jdx lr -^—e{l-x 1 )6{xi)(i)l{x l ) \x x {l-xx)\ 



dx\ 



e( Xl ) 



n-p w \ /i_ 

■ firfJ'Tr ( Tl'Pir\^ e 



dx\ 







Xi 



xl{xi 1Y «(o) 

X\ 



n-pir \ /i 







X\ 



2euv 



(130) 



For the zero-bin subtraction in the second line we kept the first term in the analog of the 
expansion in Eq. ( |123j ). Here the counterterm operator is 



o£ ] = c$ Jd V - 2 (m P ^{±m)-(ip^±}mnU;± 



;i3i) 



Pi -*o 



with 5C c t = — l/(2euv)- The tree level quark matrix element of this operator vanishes. 
Eq. ( |131| ) with 5C ct gives a vacuum to pion matrix element 



B 



ctl 



dp 2 [5{n-p^-p x -p 2 ) 4^{x u x 2 ,li 
C(0) • 



Pi ->o 



n-p n 2e uv 

This term cancels the UV divergence in Eq. ( |13(J| ) to leave the finite result 



(132) 



b w + b. 



ctl 



n-Pn { ./<< 



dxi 



Xi 



0^(0, /i) 



n-Pir Jo 



dx\ 



(133) 



As indicated the result in the first line defines the 0-distribution for this case. Again we 
should add to this the matrix element of O^} with the C^. (//_) coefficient. 

We also will need results for distributions like 0^ and (f) n but with zero-bin subtractions 
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at both ends of the integration regions 

dxdy S(l-x-y) (134) 

%0 V0 



dx r dy r — — - — 5(1— x—y) \xy\ e + c.t. 



n-p^ y ^ <%{x) qssqms | x (i_ x )| e + c .t. 
£ — gl J xil—x) 1 1 



1 

dx 







fa(x) fe(0) 
x(l— x) X 1 — x 



™-PA 2e f n x €(o)^(x-i)- 1 ^ gji^M : . t 



ik-MU (% [«o)+c(D] 

lx{l — X) X 1 — X J V |U_ / 



and 



o 



dxdy " 2 : (X ; y 2 J(l-x-y) (135) 

(a; 2 )0 (r)0 

f^) 2 ' V [dx r ^^— 2 6(x)6(l-x) |x(l-x)| e + c.t. 
V u_ / / aril — xr 1 1 

\ x ( fe) fc(0) , fal) 
\x 2 (l — x) 2 x 1—x 

v Z 1 - ' I ./l X J -co 1 — X J 

In the last line the sign for <f)' n (l) appears because we differentiate with respect to x rather 
than 1—x. The notation +c.t. indicates the matrix element of the counterterm operators 
that cancel the 1/euv divergences leaving only the logarithm. The dependence of the results 
on fi was suppressed, and terms with counterterm coefficients, (//_) [0^(0) + 0£(1)] and 
C^(n-) [0^.(0) — 0^.(1)] should be added to these amplitudes |JT . 

The final matrix elements in Eqs. ( |129| , |133|jl34|J135D have a linear ln(/x_) dependence 
which comes from the action of the 0-distribution. The coefficient of these logs is independent 
of the power taken for the x e type factors, though the analysis in section [VII B| dictates that 
the e-power should be used. The ln(/i_) dependence will be canceled order by order in a s (/z) 
by ln(/i_) dependence in the perturbative kernel J(fi-) including the coefficients C ct (/i_) of 
the counterterm operators. Just as in our scalar loop example in sections VII A 1 and |VIIBT| . 
when we consider the resulting factorization formula at the matching scale the /i + and //_ 
dependence will cancel out between logs in the coefficient functions and those in the /z± 
dependent hadronic distributions, where /i 2 = /i+yU-. 
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D. The p—TT Form Factor at Large Q 2 



In this section we consider 7*p — > tc as an example of a process with convolution integrals 
that appear to be divergent, but are tamed in our formulation of SCETn using the procedure 
in section [VII C[ The QCD dynamics are described by the p-7r form factor 



(136) 



with e 1 - the transverse rho polarization vector. At large Q 2 = —q 2 the form factor F p7T (q 2 ) 
was studied in detail in Ref. [p3H , and also was discussed in |p2| . The pion form factor 
F wn (q 2 ) and leading proton form factor F pp (q 2 ) both scale as 1/Q 2 at large Q 2 , however the 
form factor F p7T (q 2 ) scales as 1/Q 4 . This results from an additional 1/Q suppression of the 



QCD matrix element shown in Eq. ( |136| ) relative to the pion and proton cases. 

In SCET, the two LO currents that mediate any 7*Mi — > M2 transition with light 
quarks were derived in Ref. []66| . Using isospin and working in the Breit frame with incoming 
momentum transfer g M = Qn fJ '/2 — Qn fl /2, the incoming partons are n-collinear and outgoing 
partons are n-collinear. The LO matching of the QCD vector current onto SCET is 



r 



n v +n p 



doj~ 



Cl{^,UJj)Ji\p,UJj) + C 2 (p,UJ j )J^'(p,UJ j ) 



(0), 



J 



(0) 



(137) 



and involves n-collinear isodoublet Xn and n-collinear isodoublet Xn fields in SCETn. Since 



the p and 7r have opposite charge conjugation, the only relevant matrix from Ref. |R| for 
our two isotriplet mesons is 



r' 



(Qu + Qd)(r a ®T a )[ifi®i + ^75 ® ^ 7 5 ] . 



(138) 



From Eq. ( |138| ) we see that there is no 7j_, so at this order the p_i_-7r transition is forbidden. 



This is a reflection of the helicity structure of the factorization [23[ ^2, |92[ and occurs despite 
the existence of matrix elements at leading power for the relevant hadronic states: 



d n -y il±u n , x p+(p,e ± )^ = fj n-pe^S(l-x-y) (f) p± (ii,x,y) , (139) 

(n£(p') Un,uil<odn-v = ~ifn U • p 5(1 - U - v) (fi n (p,U,v). 

These twist-2 matrix elements will be useful below. The notation u n yX denotes a gauge 
invariant product of a Wilson line and quark with momentum fraction x along the quark 
arrow, u n ^ x = 5{x — 'P /n ■ p)W^n ■ Similarly, d n ,-y = 

d d) W5(y + Pyn-p) etc. Finally, the 
5(1— x— y) factors on the r.h.s. of Eq. ( |139| ) comes from momentum conservation in the matrix 
element. When implementing the 0-distribution we start with expressions containing 5- 
functions, like those shown in Eq. (|139|) . This is useful for cases where there are simultaneous 
zero-bin subtractions in variables appearing in the 5-function. 
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Once operators suppressed by 1/Q are considered, the p\_ — > n + transition is allowed p3| . 
In SCETh the relevant operators can be constructed from products of the bilinear 0(r] 3 ) 
operators from Ref. |93] 

P±(u,v) 
V&{x,y) 



1L 



A I - ^5 



d r , 



ill 



_y Pt ^ ^ P 

-1 + J I 

j -y yi y 



Tg u (u, v, w) = u n , u ^- (igB v nA _)- w d n - v , 



Vi l {x,y,z) 
A£(x,y,z) 



d n,- y ^ (wK 



d 



2 

^75 



n,-y 



(140) 



;i4i) 



where P = n a V a , V = n a V a , igB^ ± = [l/PW^in- D n ,iD^]W n ], = V£ + igB; 

and (igB^jJz = S(z — V/n ■ p)igB^ ± . As indicated, in this section we use a rescaling with 
respect to the momentum carried by the state in order to make the delta functions acting 
on the fields dimensionless. This rescaling hides the process independence of the operators, 
but makes the results simpler to present. Note that we have used slightly different notation 
for the operators than Ref. |p3| , with relations P + (u,v) = V lfw (lj), P-(u,v) = V lfw (uj), 
etc. We also included one less power of 1/V in the three-body operators. 
For the two-body operators we have matrix elements |)3] 



[tt+\P + (u,v)\o) 
[tt+\P-(u,v)\0) 



-ifnflTr S(l-U-V) </%((!, U,V) , 

— ^—^ 5{l-u-v) <f)' n a {fi,u,v) 
o 



(0\Vf + (x,y)\p+ ± (s)) = f p m p slS(l-x-y) g^(p,x,y), 
(0\V{_(x, y) |p+_(e)) = fj ^f L S(l-x-y) gft' (ji, x, y) 



2{y-x) 



8{l-x-y) g<£\p,x,y) 



(142) 



where p n = rri^ / (m u + irtd) and for later convenience we switch from the distribution g p } ' to 
a distribution g p ± . Other matrix elements are related to these by simple operator relations 
in SCET. For the three-body operators the matrix elements are |93| 
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FIG. 16: Graphs for matching onto the p-ir electromagnetic form factor. To obtain the full set 
of diagrams one must add graphs with the gluon exiting to the right (the left-right mirror images 
of c) through n) not flipping the direction of the arrows). Then to this entire set one must add 
the graphs with the current insertion on the other quark line (which can be obtained by charge 
conjugation). 



{K£(p')\Tg u (u,v,w)\0) = n-p' 



5{\ — u — v — w) 
2w 



^\Vf(x,y,z)\pt ± (P,£)) 
[0\A%(x,y,z)\p+ ± (p,e)) 



2z 



*n-pf 3o e']_ £ 



v 5(l-x-y-z) aV 



Aau\ S(l-X-y-Zj , 



(143) 



where = e MlyQ/3 n a ng/2 (which switches sign under n <-> n). In comparing to Ref. |53] note 



that we took (/ym y V) LPW = —f^fif and (f v m v A) LPW = fi p 4> 3p which agrees with the 
notation in Ref. [J23| . To compare with the other notation in Ref. ||23|| , note that (<£>^) cz = 
(^) cz = «/2, (^) CZ = 03., (^) CZ = 4> P J2, (^) cz = 9^1/2, (^) cz = 
{l-2z)g$ , {z)/4 = g£\z)/2. 

To connect the operators in Eq. (|140|) to the process ^*p\ — > 7r + we must match the full 
theory diagrams shown in Fig. [16| onto SCETn . The graphs in Fig. [I6|a-b must be expanded 
to NLO. The graphs in Fig. |l~6|c-f can be obtained from graphs a) and b) by using the tree 
level relation between the QCD and SCETn fields in place of the lowest order spinor in the 
LO part of these diagrams: 



if) = W 



1 + i^n±|] (W^ n ) . 



(144) 



Graphs g) through n) require separate computations. 

In matching these graphs onto SCETn we can make different assumptions for the scaling 
of the external lines. We work in the Breit frame and take the interpolating field for the 
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incoming p and outgoing it to be built purely out of n and n collinear fields respectively. 
The contribution that matches onto the operators in Eqs. ( |139| , |140|) will have two n-collinear 
quark fields and two n-collinear quark fields with or without an extra n or n collinear gluon. 
In addition SCETn has graphs where one or more of the above fields simultaneously become 
soft. Although we can formulate operators with soft fields in SCETn, they do not contribute 
to the p±-ir form factor in this frame. The zero-bin subtractions ensures that we will not 
double count the region of momentum space that these other operators correctly describe. 
Following Ref. [35] we note that we do not need to consider interpolating fields for hadrons 
built out of mixed soft and collinear components. These interpolating fields do not have non- 
perturbative poles as discussed earlier. Furthermore, in the Breit frame, an interpolating 
field that is purely soft would only be needed for a different physical process than the one 
we are considering (and would correspondingly require different current operators). 
To simplify the presentation we define 



x v 



S(l-x-y), 8 xyz = 5(l-x-y-z) 



(145) 



Using the computation of the tree level graphs done in Ref. [p3| , but including the ^2 p ±^ 
terms with 0-distributions as described in the previous section, gives 23 



f p AQ 2 ) 




(y 2 )0V0 



(146) 



~^~fp TTlp fir $xy&uv 
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X 0V0\ V )0 ^X0y0{U )0{V )0 ] 4 
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/J/37T 03vr(M,W,w)0 p± (x, 

y) 3uvw&xy 



2(u 2 v) (y 2 ) 2(u 2 w) (y 2 ) (uvw) y 



+ D-terms 



where 03 P = 4> 3 p A — f^ v / f^ A 4>^ p V • The "+ D-terms" factor indicates that at this order we 
must also include the renormalized coefficients Dj(/i + ,/i_) just like the coefficient C^(fiJ) 
in Eq. (pi). Here they multiply terms with p± (l), ^ } (0), ^(1), ^(0), gj£(l), 0^(0), 
(f>' w (l), and 0^(1). The /i-dependence of the distributions is suppressed for brevity. The 
range for the integrations are determined by the theta functions in the non-perturbative 
distributions, which have support from [0, 1] in their respective momentum fraction variables. 
Any variable denoted with a bar is one minus itself, x = 1 — x, etc. For the three-body 
distributions we will have two convolution integrals left after using the (^-functions, and the 



Note that we have not independently verified the calculations in Ref. 
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0-distribution must in general be treated as two-dimensional. We indicated this in Eq. (|146|) 
by having the subscript act on the product of three-body momentum fractions. The 
action of in these cases can be determined by the same steps used in section |V11 Q Due 



to the 0-distributions the result in Eq. ( J146|) is finite, independent of assumptions about the 
non-perturbative distribution functions. 



It is possible to study the /x± dependence of the result in Eq. ( |146| ). To do so one 
adopts some endpoint behavior for the distributions, and can make the action of the 0- 
distributions from the tree level jet functions explicit. A common assumption for the scaling 
behavior of the above distribution functions near their endpoints is (p w (x) ~ xx, (f) p± (u) ~ uu, 
g ( p}{x) ~ 1, gfl{x) ~ xx (so g$(x) ~ 1), <p p n {u) ~ 1, <f) 3p {x,y) ~ xyz 2 , cf)^{u,v) ~ uvw 2 . 
With this scaling behavior all integrals over three body distributions converge without zero- 
bin subtractions, and we can evaluate the two-body 0-distributions explicitly using the 
formulas worked out in section [VII Q . For the rapidity logarithms in the distributions we 
then get ln(n-p p //i_) = ln(Q//i_) and ln(n-p 7r //i + ) = ln(Q//i + ). The /i + and /i_ dependence 
in these logarithms is canceled by the /i± dependence of the _Dj(/i±) Wilson coefficients. 

In this section we showed that the SCETn zero-bin subtractions together with UV renor- 
malization yield a finite answer for the p-ir form factor at large Q 2 , given in Eq. (|146| ). 
Due to the separation in rapidity, the result has additional dependence on ln(Q) beyond 
that in the hard scattering kernel. The appearance of these logarithms is controlled by the 
powers of momentum fractions in the hard scattering kernel and the endpoint behavior of 
the distributions. At the matching scale /i + > Aq CD and ~ Q and there are no large 
logarithms. As we scale p 2 towards Aq CD keeping /i 2 = large logarithms may be 

generated. The dependence on this large log is computable, up to its normalization which 
is fixed by non-perturbative parameters. These parameters are determined once we adopt a 
model for the light-cone distribution functions. 

E. Factorization of ( B7T (E) appearing in B — > istv and B — > irir 

In this section, we study the implications of our results for B — > tt£u and the related 
process B — > titt. For small pion energy E, this process is dominated by the B* pole and can 
be studied in chiral perturbation theory . For E ^$ 1 GeV this process is also amenable 



to HQET and lattice QCD simulations [||, [97|. For large E, namely E 2 > -EAqcd > Aq CD , 
the process is expected to factor into the convolution of perturbative and non-perturbative 
contributions. 

The exploration of factorization in QCD for the process B — > i^lv has a rich history. 
In Ref. ||25|| , exclusive QCD factorization techniques were applied to the B — > tt and B — > 
uu transitions to examine the hard scattering contributions. Endpoint divergences as in 
Eq. (0) were encountered. Ref. argued that the soft endpoint contributions dominate, 



and obey heavy quark symmetry relations. This was extended to a bigger class of large 
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energy symmetry relations in Ref. |M| in the context of an EFT known as LEET ||100 



that was later shown to be inconsistent. Factorization of soft and collinear regions was 



studied in Ref. P^| , including Sudakov suppression of soft contributions, and a result for 
the hard parts was obtained in terms of leading twist B and tt distributions. B — > niu 
decays were studied with light-cone sum rules in Refs. [ |1(J1| , |102| , |103|| , and in Ref. 



a definition of the soft and hard contributions was given using the so-called local duality 
approximation. The soft part dominated numerically. In Ref. ||104| |, it was shown that terms 
with endpoint divergences could be absorbed into a soft form factor without spoiling the large 
energy symmetry relations, and a s contributions which spoil the relations were evaluated. A 
dominant non-factorizable soft form factor is an important ingredient in the power counting 
used by BBNS in discussing factorization for B — » tttt ||105|| . 24 In formulating SCET in 
Ref. @, the leading order low energy current operators for the soft B — > it form factor 



were derived, and their hard Wilson coefficients were computed at 0(a s ). Ref. ||108|| studied 
B — > ii form factors with threshold and k± resummations, argued that soft contributions 
are Sudakov suppressed in this framework, pointed out the importance of the p 2 ~ m^AQCD 
scale, and obtained a factorized result given by individual B and it wavefunctions with k^s. 
The ability to completely factorize the amplitude into individual B and tt distributions 
depending on k± is an important ingredient in the pQCD approach to B — > ixix [ |30|j . 

More recently, SCETi and SCETn have been used to analyze B — > ixiv. Refs. |35| , |109| 
pointed out that the interpolating field for the pion is purely collinear, and that the LO terms 
therefore necessarily involve time-ordered products of currents and subleading Lagrangians 
in SCETi. Thus, in the complete set of operator contributions, leading and subleading 
currents both contribute to the LO form factor. The hard-scattering and soft form factor 
contributions can be defined as subsets of these time-ordered products with and 
and it was shown that the factorization of the hard-scattering contributions involve 
a hard function for the scale ml, and also a distinct jet function for the scale mfeAQCD- 
In Refs 



75 



|, factorization in SCETn was further investigated, and it was argued that 
soft-collinear messenger contributions spoil the possibility of further factorization of the soft 
B — > 7r form factor terms. In Ref. |78|, the implications of endpoint singularities for 
factorization in SCET were investigated. It was shown that the hard-scattering contributions 
are finite to all orders in a s . It was also argued that endpoint singularities spoil possible 
factorization of the soft form factor terms below the scale mftAqcD- in Ref. ||80|| , the complete 
set of SCETn operators for the soft B — > tt form factor contribution were determined, 
including three body operators. In Refs. ||110|| a factorization formula for B — > tttt was 
investigated in SCET and it was pointed out that the same universal jet function occurs as 
in the B — > irlv process (for a precursor see ||1 1 1|| ) . 

The above results can be summarized by the following formula for the relevant B — > 71 



24 



A soft form factor does appear in the B — > Dn process [106, 107] 
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form factor, derived from a systematic expansion in E 3> Aqcd: 

U(E) = T( + \E)( b *(E)+ [ 1 dzC ( /\z,E)C^(z,E) (147) 

Jo 

rl pi poo 

= T ( - +) (E)C Bn (E) + N dz dx dk + Cj + \z, E)J(z, x, k + , E)(j)^{x)(f}%{k + ) . 

Jo Jo Jo 



Our notation follows Refs. ||109| , [L10|| , so N = /s/ 7r ms/(4i? 2 ). The hard functions and 



Cj + ^ are perturbative at ji 2 ~ m 2 , the jet function J is perturbative at jj 2 ~ -EAqcd, and 
0,r and 0^ are leading twist distributions, which are non-perturbative. Finally, ( BlT (E) = 
(m^ 2 I E 2 ')C S B ' K (E) is the unfactorized soft form factor with ( Bjr (E) containing dependence on 
the scales .EAqcd and Aq CD . At lowest order in a s (rrib) at the hard scale = 1, C { j +) = 1, 
and f + (E) = ( Bn {E) + ( Bn (E), where Cj n {E) = I dzCj v (z,E). If we also work to order 
a s (fi) at the intermediate jet scale, /i 2 ~ -EAqcd, we can expand further, and Eq. (|147|) gives 



In this result the ( Bv (E) term is left unexpanded since its factorization at scales fi 2 < EAqcd 
is (so far) unknown. In the remainder of this section we will use the zero-bin in SCETn to 
demonstrate that ( BlT (E) can be factorized further into products of twist-two and twist-three 
7r and B distribution functions. 

In SCETi the ( B% (E) term is defined by 

<7r n (p;)|T + T 3 + T 4 + T 5 + T 6 \B V ) = n-p«C(p-) C B >;) 

= 2EC{E) ( B7T {E) . (149) 

The Tj are time-ordered products of subleading Lagrangians with the leading scalar heavy- 
to-light current J® = V _ C(p-){£ n W) p - {Y^h v ). Their definitions are |TUU|, |TTU 



T = i J d'y T J(°) (0) (y) , T 3 = i J i 4 y TJ® (0) Cf q b) (y) , 

T 4 = i J*yTJ®{0) £f q a) (y) , T 5 = i 2 Jd'yd"zTJ®(0) C^(y)C^(z), 

T 6 = z 2 Jd 4 yd 4 zTJ^(0) cW{y)£$(z) , (150) 

and the presence of only J® guarantees that ( Bn (E) satisfies the symmetry relations of 
Ref. PSJ. The momentum conserving 5-function collapse the matrix element of the Tj's 
into a simple product. The momentum p~ ~ m& is large, and p~ = IE up to small power 
suppressed m n dependent terms. 

The tree level matching onto SCETn for J^ currents was carried out in Ref. [jSPR . They 
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found the position space operators 



X(0)^ x (sn) Q s {tn 



4 



Q s (tn 



few(o) 



o 



(J') 



o 



(P) 



9 



m^4cArn)x{sn)] [&(fn)^W(0)' 
X(0)^x(sn)] \Q s (tn)4 s± (un)^-H(0) 



(151) 



where the notation is defined in Ref. 



The time-ordered product T contributes to 
the matching onto operators with a gluon field strength. They also had a fifth operator, 
O5, involving a time-ordered product with soft-collinear messenger fields. In our setup the 
infrared regions associated to the messenger fields are covered by the soft and collinear fields 
so this fifth operator should not be included. With labeled fields the operators are 



4 



[LW)Jf{W%U] [(q s S)- kl ^(S%) k2 



(152) 



{Lw)Jf{w%)„ 



(q s S)^ k J^(S ] h v ) k2 
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{q s S)- k 



1 2 

ii5 



(q s S)- hl (ig^) k3 ^-(S^h v )k 2 



Tree level matching gives the Wilson coefficients |p0| 
4vra s (/i) 



Ji 



V 
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9E 2 

9E 2 
koi s (ij) 
2E 2 
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1 1 
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(153) 



uv 2 k\ ' 
w 



_(v+w) 2 vk 2 9(v+w) 2 k 2 9v(u+w)k 2 
h 8k 3 1 



+ 



_v 2 {k 1 +k 3 ) 2 k l 9v(k l + k 3 ) 2 k 1 9v 2 (k 1 + k 



In Eq. (|15|), the (igB^) was defined below Eq. QliDD , and ig& aL = [l/V S^[in-D s ,iD^ ± ]S] 
with S[n-A s ] soft Wilson lines. As in the previous section, we rescaled the delta functions 
acting on the collinear fields by n-p^ so that the subscripts involve the momentum fractions 
u,v,w, and the delta functions set either u + v = 1 or u + v + w — 1. For the soft fields 
we left the delta functions dimension-full, e.g. (q s S)_ki = (q s S)S(ki + V^). Note that the 
vacuum to pion matrix element of the operators 0% an d O2 will contribute. 25 



We thank B. Lange for pointing this out. 
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The factorization in rapidity in SCETn discussed in section [VI I| indicates that the jet 
function matching which determines Eq. ( |153| ) will give factors of a s (/x) that are perturbative 
at the jet scale fi 2 ~ £Aqcd, similar to the a s (/i) in Eq. ( |148| ). Therefore we included the 
a s (fj,) in the J[ coefficients, rather than grouping a g 2 with the operators in Eq. ( |152| ). To 
evaluate the (tt\ • • • \B) matrix elements of the 0^2,3,4, we will need Eqs. ( |139| ) and 
the pion matrix element 



,71", 



(154) 



"ft nip) 



( — u n-p) 



P + + P_)(w, v ) — i e dw T^ u iu — w,v,w 



1 \ 1 

2 lb 
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can be taken to be the decay constant in the heavy quark limit, and 



Here f B ~ rn l 

our definition for 0^ is that of Ref. | 112|| , while (j) 3 B has been studied in Ref. [|113|| and is 



proportional to — from Ref. ||1 14|| . The momentum conserving delta functions involve 
the HQET mass of the P-state, A = lim m() ^ 00 (m5 — m^). For the lowest order factorization 
formula for the soft form factor £ S7r (P), we find 
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Here v = 1 — v, the fcj-limits are —00 < k 2 < A, < k% j3 < 00 and 

8kik 2 = 5(A-ki-k 2 ), 5 kl k 2 k 3 = <$(A -ki-k 2 -k 3 ) 

S uv = 5(l-u-v), 8 UVW = 8{l-u-v-w) . 



(157) 
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The "+ D-terms" factor in Eq. (|156|) indicates that at this order we must also include the 
renormalized coefficients _Dj(/i + ,/i_) for operators like the ones in Eq. ( |128| ). In Eq. ( |156| ), 
we have obtained a finite factorization formula for the soft B — > 7r form factor. The structure 
of this contribution is very similar to the hard-scattering term in Eq. (|148 ), except for the 
necessity of the zero-bin subtractions. 

To study the /x± dependence in Eq. ( |156| ), we can adopt the standard endpoint behavior 
for the distributions. For simplicity, we can adopt a behavior for 03b(/ci, k 2 , k 3 ) so that its 
tree-level integrals also all converge without applying the zero-bin and renormalization. This 
leaves ^{u) ~ uu, <j)%{k\) ~ k%, 4>s(ki) ~ 1, and (p 37T (u,v,w) ~ uvw 2 . The zero-bin is then 
required for cj) n (u, v) / (v 2 )0 given in Eq. ( |129| ) and for (j)^{u, v, w)/[v 2 w(u+w)]0 which gener- 
ate ln(/x_/p^) = ln(/i_/2i?) terms, and for k 2 )/(kl)0 and 4>g(ki, k 2 )/{ki)0. Following 
the usual procedure, we find that these last two terms each generate a ln[/i + /(n-i>A)] 

At higher orders in a s , we anticipate that the factorization formula will take the form 

C Bn = N Jdu Jdki Ji(u, k i ,E,n,ij ± )(f) n (u,n,fi-)(f) B (k 1 ,n,n + ) (158) 
+ N Jdu Jdk x J 2 (u,k 1 ,E,fi,^ ± )((j)l + ^(j)' 7T a )(u,fi,fi-)(f)^( y k 1 ,fi,fi + ) 
+ N Jdudv y"d/ci J 3 (M,w,fci,E,/i,/i ± )0 3w (u,w,/i,/i_)0+(A; 1 , y u,/i + ) 
+ N du dhdk 2 J4(u J k ll k 2j E J ij, J fi ± )(f) 7T (u J iJ,fi-)(f) 3 B(kuk2 J n,iJ, + ). 



where we have used the etc. functions to reduce the number of integrations, and we 
note that the zero-bin subtractions are present on the remaining variables. Here the fi± 
dependence in the <j>^s and 0^'s is a short hand for the ln(/z_) and ln(/x + ) terms that are 
generated from the 0-distribution depending on the structure of the Jj's and the endpoint 
behavior of the distributions. The completeness of the mixed soft-collinear basis of operators 



Oi,2,3,4 found in Ref. [§0J should guarantee that it is the non-perturbative functions shown 
in Eq. ( |158| ) which will show up at any order in a s in the matching. 26 One-loop corrections 
to T^ + \ Cj, and to J are known 0, [H], |115| , |116|| . One-loop corrections for the Jj can be 



computed following the method outlined here. 

Since the factorization formula for B — > ttti involves the same ( Bn and n functions, 

A(B -> vrvr) = X^A^ + ^^{f^^ j} u [T lc (u)+T 2C ( M )j 0» 

+ U Jdu Jdz [T 1 j(u,z) + T 2J {u,z)] ( f>*(u)(f*{z)y (159) 

the results in this section immediately carry over to that process. In Eq. (|159|) = 
Cfj)(E = m b /2) and we quoted the result from Ref. [|110|| . Our result in Eq. (|156|) lends 
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The results in Eq. (77,78) of Ref. M, allow for the elimination of 4>p and 4>'^ in terms of fo^ if so desired. 
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support to the power counting for B — > tttt used in Ref. [ 11011 , where ( Bn and (f n were 
treated as being parametrically and numerically similar in size. It also indicates that the 
entire non-leptonic tree amplitude for B — ► ttti can likely be written in terms of individual 



B and 7r distribution functions. This differs from the BBNS [105, TT7\ type factorization 



where the soft-form factor is taken to be non-perturbative. Our results appear to indicate 
that the soft-form factor factorizes, but with the SCETn factorization in rapidity space that 
is different from the standard type of factorization formula. A factorization of the form 
factor in terms of individual B and tt objects is similar to the result found in the pQCD 
analysis in Ref. [^, [30|, however we do not find kj_ dependent functions. A more detailed 
study of Eqs. (|156|) will be reported on in a future publication. 

VIII. CONCLUSION 

This paper discusses how to formulate effective field theories with multiple fields for the 
same physical particle, where each field represents a region of momentum space with a 
different power counting. This allowed for a tiling of all the infrared momentum regions 
of the theory. The technique was applied to a non-relativistic field theory (NRQCD) and 
to field theories for energetic particles, including SCETj with collinear and ultrasoft quarks 
and gluons, and SCETn with soft and collinear quarks and gluons. Effective theories with 
multiple fields for the same particle can be formulated using fields with labels that distinguish 
the different momentum modes. In converting the sum over momentum labels into an 
integral, one must be careful with the zero-bin, where the momentum label is zero, since 
this corresponds to a different degree of freedom. Here the proper treatment of the zero-bin 
was investigated in detail. Sums which do not include the zero-bin are converted to integrals 
over all space with a zero-bin subtraction. The zero-bin subtraction removes the support 
of the integrand in the infrared region of momentum space where overlap between different 
modes could occur, and thus avoids double- counting between different modes in the effective 
theory. 

The zero-bin subtractions give a definition of NRQCD and SCET independent of the 
choice of UV and IR regulators, allowing the use of regulators other than dimensional reg- 
ularization if desired. They also solve a number of puzzles encountered in the literature on 
collinear factorization in QCD and non-relativistic field theories, which are associated with 
unphysical pinch and endpoint singularities, as well as puzzles in the NRQCD and SCET 
literature with UV and IR divergences that occur at intermediate length scales. 

In NRQCD, the zero-bin subtractions eliminate pinch singularities in box-type graphs by 
properly distinguishing between potential and soft heavy fermions. This result applies to any 
non-relativistic effective theory (for example, it also simplifies the resolution of the puzzle 



discussed in Ref. [|118|| involving soft pions interacting with nucleons). Zero-bin subtractions 



also distinguish between soft and ultrasoft gluons in NRQCD. In graphs with soft loops, 
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they convert infrared divergences into ultraviolet divergences. Converting an infrared into 
an ultraviolet divergences in the soft sector corresponds to a rearrangement of degrees of 
freedom that was known to be necessary in NRQCD, and had been implemented by hand 
as the "pull-up" procedure |TJ|, [I5J. Here a proper treatment of the zero-bin enabled us to 
derive the pull-up mechanism directly from the effective Lagrangian. 

In SCETi zero-bin subtractions occur in collinear loops and avoid the overlap with the 
ultrasoft momentum region. These subtractions remove infrared divergences and induce 
new divergences in the ultraviolet. Previous results in the SCET literature used anomalous 
dimensions computed by including infrared divergences in collinear graphs as though they 
were ultraviolet divergences. This gives (by hand) the correct result for the anomalous 
dimensions. The zero-bin subtractions allows the computation of anomalous dimensions in 
the effective theory to be carried out by the standard renormalization procedure, in terms of 
the counterterms used to cancel ultraviolet divergences, and also ensures that IR divergences 
of the full theory are properly reproduced independent of the regulator choice. In inclusive 
decays at large energy, the zero-bin subtraction applies to both virtual loop integrals as well 
as phase space integrals for real emission as described in section |V1 D[ We believe that zero- 
bin subtractions will play a role in the SCET matching calculations for parton showering 
carried out in Ref . ||119|| , and perhaps also to the subtractions carried out in Ref . ||120| . They 



are likely to have implications for the singularities encountered with /^-distributions with 



light-like Wilson lines which were discussed in Refs. |3l], [32| . 

In exclusive decays, the zero-bin subtraction gives factorized decay rates and cross- 
sections with finite convolution integrals. Convolution integrals of the perturbatively calcu- 
lable kernel with hadron wavefunctions are sometimes naively divergent at the endpoints, as 
is the case for the p-7r form factor, and in factorization of the soft form factor in B — > niu 
decays. The zero-bin subtraction avoids double counting between soft and collinear modes 
in SCETn, and this converts the unphysical infrared divergence in convolution integrals into 
an ultraviolet divergence. These ultraviolet divergences are canceled by operator renormal- 
ization. The final convolution integral is finite, determined by a distribution we called 0. 
The 0-distribution is a plus-type distribution augmented by additional non-analytic hi(E) 
dependence induced by the renormalization. Independent of the choice for the UV and 
IR regulators, our formulation of SCETn has only soft and collinear degrees of freedom 
and does not have soft-collinear messenger modes [75J as explained in section |V11 A 1| and 
Appendix [B]. 27 

Thus one gets finite unambiguous formulas for amplitudes in SCETn, free of endpoint 
singularities in the convolution integrals. In cases where there would have been a divergent 
convolution we get a separation of modes in rapidity space with variables fi±. This preserves 



27 In Refs. KQ, u\ 121 1 SCET was reformulated in position space to avoid the sums over label momenta. To 



formulate the zero-bin subtractions in position space one can take the Fourier transform of Eq. ( |l7|) . 
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the power counting and naive factorization in the kinematic variables that would be present if 
the convolutions were finite. To illustrate this, we derived finite amplitudes for the p-ir form 
factor at large Q 2 , and for the soft form factor function ( B7T (E n ) that appears in B — > tc£u 
and B — > tttt decays. Without a proper treatment of the zero-bin the convolutions in these 
factorization formulas would suffer from endpoint singularities. Applying the technique to 
the soft form factor for B — > tttt decays allowed us to derive an amplitude for non-leptonic 
decays that is entirely in terms of B and 7r distribution functions at lowest order. 

In previous literature, the field theories, SCETi and SCETn, have provided an algebraic 
means of deriving factorization formulas at any order in Aqcd/Q, but not the means to 
guarantee that the manipulations would result in finite factorization formulas. Our work 
suggests that this will indeed be the case, so that the power expansion of observables can be 
carried out to any desired order without encountering singularities in convolution integrals. 
It still remains to carry out the full derivation of an SCETn factorization formula and explore 
the renormalization group properties of the resulting amplitudes - tasks which we leave for 
future work. The zero-bin procedure provides the freedom to tile the infrared of an EFT 
with suitable degrees of freedom, and makes the connection between the choice of degrees of 
freedom and the power counting expansion clear. Exploiting this, we converted the question 
of finding a complete set of degrees of freedom for an EFT at any order in the expansion, 
to a question that is easier to answer physically, that of identifying the relevant operators 
that occur at leading order for an observable, and give a proper formulation of the power 
counting. 
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APPENDIX A: RAPIDITY CUTOFF LOOP INTEGRALS IN SCETn 



In this section we give some details on the calculation of the integrals required for the SCETn diagrams 
in section VILA. This calculations uses dimensional regularization for IR divergences, and cutoffs on a Wick 
rotated rapidity variable, £L to regulate rapidity effects in the UV as in Eq. (p4h. The integrals we wish to 
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compute are 
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with £+ > and p- > 0. We use variables {k+Xk} for the soft integral with Wick rotation = i£' k , and 
{k~,(k} for the collinear integral with £/. = — «Cfcj as discussed in section VII A It is easy to verify for I s 
and I c that these Wick rotations about the origin do not encounter any poles. 

Soft Integral in SCET n 



Consider the soft integral in Eq. (|A1|) and perform the k± integral: 
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Let A; = ^k+ , and note that dk = \k+\d(^ once we integrate — oo < k+ < oo and — oo < C < oo. Thus 
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Rescale by a positive constant k+ = (£+x)k'+ and for simplicity rename k'+ = k+. Then rotate £ = i£', so 
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We divide the k+ integral into k+ > 1, < k+ < 1, and k+ < 0. This gives 
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In the second equality note that the change of variables k+ — > 1 — fc + for fc + < causes the second and 
third terms to cancel. Putting the pieces back together and multiplying by the /i 2e this gives 
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Now set Cmin = an d Cmax = ° 2 an d note that (±i) e = exp(=Fz-/re/2) , so that we have 
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(In the last line when expanding in e we have multiplied by the exp(e7E) factor to put fx into the MS scheme, 
while the necessary (47r)~ e factor was removed already in Eq. (A2).) Eq.(A7) is the result quoted for the 
soft computation in Eq. (p8|) in the text. 



Collinear integral in SCETn 



Lets repeat the computation in the previous section for the collinear integral in Eq. (97). The integral is 
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Rcscale by a positive constant k + = (p x)k' + (for simplicity renaming k' + — > fc + ), and then rotate £ = — 
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where in the last line we noted that the /c integral is identical J(l/£') defined via Eq. (A5). Multiplying 
by /z 2e this gives 
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Now we add the contributions from the regions (' 6 [a 2 , oo] and £' € [—oo, —a 2 ] to give 
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(A12) 



This is the collinear integral result quoted in Eq. (|99|). A simple way to get this answer is to note that the 
original collinear integral is identical to the soft integral with the replacements k + «-> k~, p~ «-> and 
a — > 1/a. The answers in Eqs. (A7) and (A12) agree with this. 
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APPENDIX B: SCET n LOOPS IN DIM. REG. WITH DIFFERENT IR REGU- 
LATORS 



In this appendix we repeat the matching computation done in section VII B 1 of a scalar loop integral in 
SCETn- We use dimensional regularization, but modify the treatment of the IR regulator. The structure of 
the full theory and effective theory diagrams changes, but again the IR divergences are properly reproduced 
and the same contribution to the matching coefficient is obtained. The calculation is done for two classes 
of IR regulators: i) taking three different IR masses, mi, m%, and 7773 rather than just the single mass used 
in section VIIBI, and ii) with mi ^ 0, mi = 777,3 = 0, and external momenta offshell, £ 2 7^ and p 2 7^ 0. 



Finally, in a part iii) we discuss subtleties related to the mi — limit of these two cases. For simplicity we 
leave off the diagram prefactor ieg 2 G/(p~ £ + ) and just quote results for the integrals in this appendix. In 
all cases p~ > and £ + > 0, and (p~£ + ) is the perturbative scale. 



i) Three IR Masses, mi, ni2, HI3 



The full theory loop is the generalization of Eq. ( VII B 1 ) with three IR masses, 
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In Eq. (Bl) factors of the IR regulators, m\, 777J, and m\ appear in all propagators, external momenta are 
taken onshell, p 2 = I 2 — 0, and we have expanded in mf / (p~ £ + ) . The result is valid as long as {p~£ + ) ^> m 2 , 
(p~£ + )m 2 3> m\m\, and can not be used for the case mi = since it blows up. The result which is valid 



for 771! — > and also reproduces Eq. (Bl) is given below in Eq. (B22). The mi — result is in Eq. (B24) 



The LO SCETn currents for dimensional regularization are given in Eq. (105). Using the 777,1,2,3 IR 
regulators for the scalar and collinear loops in Figs. [l5|b,c we have 
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Note that we keep m| and m\ only in the propagators that are allowed to become small by the power counting 
in SCETn ■ Eq. ( |l7| ) tells us that we have zero-bin subtractions for the soft and collinear diagrams which 
avoid the IR singularities from the [— p~k + ] and [— £ + k~] propagators. The naive integrals and subtraction 
integrals are 
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To determine the form of the subtraction integrals, we considered the collinear limit of the soft loop mo- 
mentum in /loft, and the soft limit of the collinear loop momentum in /^ alar . Note that the m 2 and m 2 
dependence is dropped in the subtraction integrals because in these limits £ + k~ ^> mf. and p~k + ^> m 2 . The 
UV rapidity regulator factors, | • • • | e are not affected by the subtractions (despite the way we are organizing 
the computation here, one really has an integrand defined with subtractions and then multiplies it by these 
factors) . From Eq. (|l7|) the differences I^ft™ - ^osoft 1 an d ^cn alar _ ^ocn' ar wm g ive tne result for I^ av and 
jscaiar respectively. 

For the soft graph we do the k~ integral by contours. Due to the pole structure this restricts the k + - 
integration to the region < k + < £ + . The k± integral is then done. For the soft subtraction integral we 
follow the same procedure which this time leaves the integration region < k + < oo. We find 
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The double counting with the collinear integral comes from the k + — > part of the integral, but the 
divergence from this limit exactly cancels in I^f^ 1 = -f|oft ar — ^OsoftT as l° n § as 1711 ^ ®- Computing the 
integrals we find 
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For the collinear integrals we do the contour integration in k + which restricts the remaining integration 
region in k" . For the naive and subtraction integrals we find 
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The subtraction integral cancels the singularity in I[ 



calar 



as k 



as long as mi ^ 0. The complete 



collinear result, J^ alar = J^ alar — JgJ: alar , is very similar to the soft result 
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The results in Eqs. (B5) and (B7) have 1/euv hi(m 2 ) divergences, which are canceled by the <fi(0, //)/e type 



counterterms. For this scalar calculation these divergences are canceled by a graph containing the insertion 
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of the renormalized currents in Eq. (112) with additional counterterm coefficients for the convolution integral 



as given in Eq. (|113| ). Contracting the scalar gluon as in Fig. |15p), using the same IR mass regulator, and 
pulling out the same prefactor as the other diagrams gives 
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Note that dependence on 777,2. 3 drops out of the answer for I^ alar . Due to our choice of SC ( - od ' ^ this exactly 
cancels the 1/euv m (TO 2 ) terms in the collinear and soft loops. Adding the soft, collinear, and counterterm 
graphs we find the full SCETn result 
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The effective theory result in Eq. (B9) has UV divergences which are the same as in Eq. (115), and are 
canceled by a counterterm for the jet function coefficient j( 0a \ as given in Eq. (116). The renormalized 
EFT result is 
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The first three terms exactly reproduces the IR divergences in the full theory result in Eq. (Bl), including 
the entire functional dependence on the ratios of to 2 , and the fourth term vanishes since /i 2 = /i + /z_ . The 
difference of the remaining finite terms gives a contribution to the one-loop matching coefficient 
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This result exactly reproduces the matching coefficient in Eq. (119), as anticipated. In the limit 7712,3 ~~ * 



all results go smoothly over to those in section VTIB 1 



ii) Offshellness p 2 = -P 2 ^ 0, £ 2 = — L 2 ^ 0, with mi/0 and m 2 ,3 
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The full theory loop integral is now 
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where £ 2 = —L 2 , p 2 = —P 2 and we have expanded in P 2 /(p~£ + ), L 2 /(p~£ + ), and m 2 /(p~£ + ). The result 
is valid as long as (p~£ + )ml 3> P 2 L 2 , and so can not be used for the special case toi = 0. The result for 
toi — > is discussed below in case iii). 
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Using the same IR regulators for the scalar and collinear SCETn loops in Figs. |15|b,c we have 
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Note that we keep I 2 and p 2 only in the propagators that are allowed to become small by the power counting 
in SCETn ■ Eq. (|l7|) tells us that we have zero-bin subtractions for the soft and collinear diagrams and the 
naive integrals and subtraction integrals are 
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To determine the form of the subtraction integrals, we considered the collinear limit of the soft loop mo- 
mentum in 7|oft lar , and the soft limit of the collinear loop momentum in i^ I c 1 alar . The p 2 and £ 2 dependence 
is dropped in the subtraction integrals because in these limits £ + k~ 3> — £ 2 and p~k + 3> — p 2 . Note that 
the subtraction integrals are identical to the case with mi, 2, 3 7^ 0. The final results I^of^ and /scalar are 
defined by the differences / s S o f alar - ^sofT an d C' ar ~ ^ocn a ' respectively. 

To compute the soft graph we work in the frame where I = 0, thus £ + £~~ = —L 2 < with £ + > 0, so 
the offshell momentum is l~ = —L~ < 0. We do the k~ integral by contours. Due to the pole structure 
this restricts the A: + -integration to the region < k + < £ + . The k± integral is then done. For the soft 
subtraction integral we follow the same procedure which gives the integration region < k + < 00. Thus 
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The double counting with the collinear integral comes from the k + 
divergence from this limit exactly cancels in /|^ lar 
integrals we find 
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— ► part of the integral, but the 
as long as mi ^ 0. Computing the 
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For the collinear integrals we take p- 1 = 0, so p p + = —P 2 < with p > 0, and it is p + = —P + < that 
takes the p 2 offshell. We do the contour integration in k + which restricts the remaining integration region 
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in k . For the naive and subtraction integrals we find 
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The complete collinear result, /^ alar = 4 s S alar — ^ocn lar > i s ver y similar to the soft result 
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The results in Eqs. (B16) and (B18) have l/euv m ( m2 ) divergences, which are canceled by the <p(0,fi)/e 
type counterterms just as in our 771^2,3 case. The countcrterms and result to be added are the same as 
in Eq. (B8) and exactly cancels the l/euvhi(m 2 ) terms in the collinear and soft loops. Adding the soft, 
collinear, and counterterm graphs and simplifying we find the full SCETn result 
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The effective theory result in Eq. (B19) has UV divergences which are canceled by the counterterm for the 
jet function coefficient J^ a \ as already given in Eq. (116). The renormalized EFT result is 
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The first three terms exactly reproduces the IR divergences in the full theory result in Eq. (B12), including 
the entire functional dependence on P 2 /m 2 and L 2 /m\. The fourth term vanishes since [i 2 = /i+ The 
difference of the remaining finite terms gives a contribution to the one-loop matching coefficient 
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This result exactly reproduces the matching coefficient in Eqs. (HE) and (Bll) using a different IR regulator. 
This is as expected since the full theory was UV finite, and the same UV regulator was used in the SCETn 



calculation. In the limit L 2 , P 2 — > all the results go smoothly over to the results presented in section VII B 1 



iii) The limit mi — > of cases i) and ii) 

Finally, we discuss the limit mi — > of the IR regulators considered above in cases i) and ii). This is 
not a smooth limit in either the full or effective theories. In the following we use the notation Q 2 = (p~f + ) 
as a shorthand for our large perturbative scale. 
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We first consider the full theory loop integrals, but in expanding out the IR regulators we keep the first 
subleading terms in the expansions in cases where the leading term vanishes as m 2 — > 0. So for the expansion 
in m 2 2i3 /Q 2 in case i), we keep subleading m| 3 terms if the leading term is proportional to m\. For the 
expansion in m 2 /Q 2 , L 2 /Q 2 , and P 2 /Q 2 in case ii) we keep subleading L 2 and P 2 terms when the leading 
term is proportional to m\ . For the TOi,2,3 7^ regulator this gives 
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while for the m 2 7^ 0, p 2 
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From Eq. (B22) we see that as long as Q 2 m 2 m 2 m§, as is the case if all the IR masses are the 



same order in the power counting, then £ — > Q m 2 , and expanding Eq. (B22) reproduces the result quoted 
If we set mi = 0, then the leading m 2 /Q 2 terms in the double log and di-logs vanish, and 



in Eq. (Bl 



subleading terms regulate the IR divergences. In this case £ 



oTnf;, and we obtain 
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We have checked that Eq. (B24) agrees with the result obtained by setting mi = before evaluating 
the integral. Since subleading terms are regulating the IR divergences the result depends on the product 
ln(m|)ln(m|) and thus no longer has a form that can be factorized in a straightforward manner into soft 



and collinear parts. The situation is very similar for case ii). Expanding Eq. (B23) when Q 2 m\ ^> L 2 P 2 
gives £ — ► Q 2 m\ and reproduces the result quoted in Eq. ( B12 ). If we set mi = then the leading term 
in £ vanishes and subleading L 2 and P 2 terms regulate divergences in the double log and di-logs, with 
(, —> —L 2 P 2 . In this case 
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To see why this happens we can examine the IR divergences in the full theory integral. Divergences 
occur for k — > p~, k — > t + , and k — > 0, and the issue with mi — > arises from the k — > case where 
the propagators carrying the soft and collinear momenta are both singular in opposite light-like directions. 
limit of the denominator of the integrand in Eq. (Bl) is [— l + k~ — m 2 ][fc 2 
L 2 ][k 2 



The k 



m 2 ][— p k^ 



and for Eq. ( |B12j ) is [— £ + k~ — L 2 ][k 2 — m 2 ][— p~k + — P 2 ]. The effect of the IR regulators is pictured by 
the solid curves in Fig. [1^, where we show how they shield the integrand from blowing up when we go 
towards the k~ = or k + = lines. Without the mi regulator the intersection of the k + = m 2 /p~ and 
k~ = m\ji + lines generates unphysical sensitivity to a very small scale ~ m\m\j Q 2 . For the case with 
offshellness it is the intersection of the lines k + = P 2 /p~ and k~ = L 2 /£ + generating sensitivity to the very 
small scale ~ L 2 P 2 /Q 2 . The sensitivity to this new small scale was first pointed out in Ref. |Q where 
the IR regulator L 2 7^ 0, P 2 ^ was used. The result in Eq. ( B25 ) agrees with the one studied there. In 
Ref. [ [75] messenger modes with very small invariant mass were added to the effective theory to account for 
the dependence on this "messenger scale". However, in QCD the sensitivity to this small messenger scale 
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FIG. 17: Regulation of IR divergences in SCETn in perturbation theory in the k + -k~ plane. 
The c n and s denote the collinear and soft modes respectively, and the dashed k + /k~ = /7 + //i~ 
line indicates how these modes are distinguished in perturbation theory. For k^ — > the solid 
lines indicate at what scales the IR divergences are cutoff by the m 2 , mf,, and 777.3 regulators (and 
also for the mf, L 2 , P 2 regulator choice). Since p~ S> 1x1% the intersection of the k + = m 2 /p~ 
and k~ = m 2 ./^ lines is always below the curve k + k~ = m 2 (and same for k + = P 2 jp~ and 
k~ = L 2 /£ + ). This intersection occurs at the messenger scale. 



is unphysical because IR divergences are cutoff at an earlier stage by Aqcd j85|. In perturbation theory 
with mi 7^ 0, sensitivity to the messenger scale also never appears. Comparing Fig. [l3| with Fig. [l7| we see 
that the m\ ^ regulator behaves in a similar manner to Aq CD . Other IR regulators are also known which 
remove the unphysical sensitivity to the messenger scale, including analytic regulators Wm, and an energy 



dependent gluon mass |65 



It is possible to choose IR regulators that complicate the choice of the matching coefficient; m\ = is 
such a choice. We have seen the complication in the full theory result in Eq. ( B22 ). In the effective theory 
the choice mi = also causes problems for I^ al — ^osoftT an d 7j?£ alar — /ocn' ar - in this situation the k± 
integration for the subtraction integral is scaleless, and the subtraction terms do not cancel the problematic 
k + — > and — > regions in the naive integrals. For example the naive and subtraction integrals in 



Eq. (B4) scale with a different power of k + as k + — > 0. The same is true in for Eq. ( |B15| ), and for the 
collinear integrals. If one chooses mi = some of the IR divergences are regulated at the scale m\m\/Q 2 or 
L 2 P 2 /Q 2 , and so part of the IR regulator is of subleading order in the power counting with our definitions of 
the soft and collinear modes. This makes it difficult to compute the IR behavior in the effective theory, since 
one has to sum up a class of subleading terms in the power counting to all orders. We expect a resummation 
procedure in SCET n could be developed to reproduce the matching result in Eq. (11£) with mi = for the 
IR regulators considered in i) and ii), but we will not attempt it here. 



It is worth emphasizing that the matching coefficient, Eqs. ( |119| Bll ), which depends on the difference 
between the full and effective theory results, does not depend on the IR regulator, as we have demonstrated 
with several different calculations. The matching coefficient is only sensitive to ultraviolet physics. It can 
be computed by using an IR regulator that is homogeneous in the power counting. We have seen this 
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explicitly by computing the matching for arbitrary 7711,2,3, and for P 2 1 L 2 ) m\ with mi 7^ 0. The purpose 
of the perturbation theory Feynman graph computations in the full and effective theories is to compute 
the matching coefficient Eq. (B4), which is a short distance quantity. The effective theory is not being used 
to reproduce the IR behavior of the full theory in perturbation theory for arbitrary IR regulators, especially 
regulators that only become effective at a small scale that is subleading order in the power counting. The 
effective theory is being applied to QCD, where the infrared dynamics is nonperturbative, and cutoff at 
Aqcd- 
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